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ABSTRACT 


Based on the experimental observation that in the fully 
developed flow region of a submerged three dimensional jet, the 
velocity profiles in the normal and in transverse planes are all 
similar, a suitable form of nondimensional velocity distribution was 
introduced. The boundary layer type of equations for the three di- 
mensional free jet was then analysed and relationships for the maxi- 
mum velocity decay and growth of length scales were established. The 
approach was then extended for a three dimensional wall jet growing 
on a smooth wall in a very wide channel, assuming that the forward 
momentum of the jet was preserved. Flow from deeply submerged out- 
lets of various geometric shapes resting on the bed, to discharge 
tangentially were investigated experimentally. For each shape of 
outlet, the velocity profile normal to the wall as well as in the 
transverse planes were found to be similar. The centerline non- 
dimensional profile normal to the wall was found to be described 
closely by the classical wall jet curve, whereas, the transverse pro- 
files were described by Goertler's curve for a free circular jet. 

The streamwise variation of the velocity and length scales were 
studied and compared with theoretical laws. From dimensional analysis, 
the square root of the area of cross-section of the outlet was found 


to be a characteristic length parameter and was found to be useful in 
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correlating the experimental data for maximum velocity decay, regard- 
less of the outlet shape. The transverse distribution of bed shear 
was investigated and was found to be similar. Curves for estimating 
the centerline bed shear are presented in convenient forms. 

Other aspects of practical interest, such as, the effect of 
the side walls of the channel and the bed roughness, on the growth of 
a circular wall jet were also investigated. A circular wall jet with 
Swirl was then studied experimentally and was found to be a very 
effective means of reducing the maximum forward velocity in comparatively 
shorter distances. Further, the similarity of velocity profiles and 
the bed shear distribution for the swirling jets were investigated. 

The similarity analysis for velocity distribution was found 
to be useful even in more general types of flows. The turbulent 
compound annular shear was analysed to establish a relationship for 
the growth of the half velocity length scale. Experiments were con- 
ducted for four different secondary stream to jet velocity ratios, 
and the results justify similarity of velocity profiles and the growth 
law. The constant of proportionality in Prandt1l's mixing length re- 
lationship was found to increase with the velocity ratio, and was far 
from being a constant as assumed by other investigators (Appendix A). 

In Appendix B, some of the available data on three dimensional 


free jets were analysed for velocity decay and growth rates. 
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CHAPTER I 
INTRODUCTION 


1.1 General 
i In a number of hydraulic structures a high velocity stream 
meets an ambient fluid at rest and gets diffused. As the high velocity 
Stream meets the stationary water, the so called surfaces of separation 
or discontinuity would appear. Tangential separation surfaces are 
highly unstable and result in eddies that move in a random fashion 
both along and across the stream, thus bringing about an exchange of 
momentum between adjoining layers. This results in a region of finite 
thickness having a continuous distribution of velocity and is termed 
a shear layer. In most practical cases the Reynolds number of flow 
is large and the shear layer is turbulent. This turbulent shear 
layer gradually extends into the non-viscous core of the oncoming jet 
of fluid and beyond a comparatively short distance the shear layer 
extends over the entire thickness of flow. 

The type of turbulent jet, that has been studied most, is 
the one diffusing through a medium at rest. The two common types 
investigated in detail are the two dimensional or the plane jet, and 
the axially symmetric jet. If the jet diffuses in the absence of 
any solid boundary, it is called a free jet and if it diffuses near 


a solid boundary, it is known as a wall jet. A frequently encountered 
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2 
example of a wall jet is the one with a jet of fluid impinging tangentially 


on a solid wall and then growing on it. 


1.2 Hydraulic Outlets 


Control gate structures are frequently encountered in any 
irrigation and power development project. During certain periods, these 
control structures operate under submerged tail water conditions and 
in general the flow from any submerged outlet can be analysed as a 
turbulent jet. If the boundaries are sufficiently away, the flow 
can be treated as a free jet problem and the information available 
on it is quite extensive (Abramovich, 1963). However, in practice, the 
flow from the outlet issues close to a boundary, normally the bed of 
a channel, and has to be analysed as a wall jet. 

In few cases the outlet might occupy the full width of a 
channel and could then be treated as a plane wall jet. Considerable 
amount of information is available on plane wall jets (Schwarz and 
Cosart, 1961; Meyer et al., 1961) and it has been demonstrated 
(Rajaratnam and Subramanya, 1967) that the plane wall jet serves as 
an useful model to describe the flow below sluice gates in hydraulic 
engineering. However, in most cases the flow emanates from outlets 
that do not occupy the entire width of the channel and have finite 
aspect ratios; being situated right on the bed of the channel or at 
a small height from the bed. Some information on diffusion of wall 
jets into wider channels is available (Rajaratnam and Subramanya, 


1967), but the investigation does not include any shape other than a 
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3 
rectangle. A preliminary experimental study is available on flow from 


non-rectangular outlets (Rajaratnam and Muralidhar, 1969). 


1.3 Existing Work 


The experimental investigations on jets of finite aspect 
ratio or the three dimensional jets are limited. Yevdjevich (1966) 
studied the diffusion of three dimensional free jets issuing from 
rectangular slots having a wide range of aspect ratios. Sforza et al. 
(1966) and Trentacoste and Sforza (1967) studied the behaviour of 
free jets emanating from rectangular nozzles having different aspect 
ratios as well as from other nozzle shapes that include a triangle and 
an ellipse. The investigation on a three dimensional wall jet was 
more restrictive and covered only the rectangular nozzle (Viets and 


Sforza, 1966; Sforza and Herbst, 1967) having a few aspect ratios. 


1.4 Present Investigation 


An understanding of the flow field and shear stress distribu- 
tion on the channel bed resulting from the operation of submerged 
culverts and gated structures will be of great practical use, particu- 
larly in designing efficient energy dissipating devices and bed pro- 
tections. Information on the behaviour of diffusion of wall jets 
issuing from triangular, elliptical and the more frequently encountered 
circular shape outlets is not available at present. In addition the 
available literature on rectangular wall jets is very limited and many 


aspects of the problem remain unanswered. The present investigation 
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was undertaken in an effort to understand the mean flow characteristics 
for a variety of shapes of practical interest and develop relation- 
ships that would apply to an outlet irrespective of its shape. Further- 
more, the experimental undertaking is also an effort to accumulate 
sufficient information, which would enable to predict the flow charac- 
teristics for engineering design purposes with some degree of confidence. 

In the experimental program, to start with, a simple case of 
a three dimensional wall jet diffusing in a wide channel and a smooth 
wall with no swirl was studied. However, this ideal case is far from 
reality and in most practical cases the channel side walls are close 
to one another. Proximity of the side walls would interfere with 
the growth of the jet and might significantly alter the behaviour of 
the jet. In an attempt to study this behaviour, the experimental 
program was extended to investigate the diffusion of circular wall 
jets in channels of finite width. 

Another significant aspect in which the ideal case considered 
so far differs from practice is the roughness of the bed or wall. 
In reality, under most hydraulic outlets the bed is rough. To assess 
the order to which the ideal case is getting affected by the rough- 
ness elements, the present investigation includes the case of a 
circular outlet on a rough bed. 

In a free circular jet, a certain amount of swirl or rotation 
before it comes out of the nozzle helps the jet to spread more rapidly, 


and the velocity field decays much faster (Chigier and Chervinsky, 1965). 
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It is therefore thought that giving a rotation to the circular jet 
might greatly assist in diffusion of the jet efficiently in a shorter 
reach and would be useful in reducing the length of apron or protection 
needed below the outlet. To explore this aspect the investigation 


was further extended to diffusion of circular wall jets with swirl. 


1.5 Presentation of the Program 


The following is an arrangement adapted in presenting the 
details of the present investigation. 

A theoretical consideration is presented in Chapter II. 
Based on the observations that the velocity profiles are similar, the 
theoretical analysis predicts certain basic relationships for the 
decay of velocity and the rate of spreading of the jet. Dimensional 
considerations are also presented to arrive at a number of useful 
dimensionless parameters, which set the stage for presentation of the 
experimental results. 

Chapter III gives the details regarding the experimental 
arrangements and the measuring techniques adapted. Certain special 
types of probes had to be built and tested to measure the three 
dimensional velocity field and the boundary shear distributions. The 
calibration techniques and the methods of using the calibration charts 
are also discussed. The limitations of the measurements are then 
described. 

In Chapter IV, the investigations on three dimensional 


turbulent wall jets diffusing in a wide channel and on a smooth bed 
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with no swirl arepresented. The shapes studied include a circle, a 
Square, anequilateral triangle, an ellipse and rectangles of different 
aspect ratios. Tne results are presented in convenient forms that 
would predict the velocity decay, the growth of the jet and the bed 
shear distribution in the channel. 

Investigations on the effects of side walls on circular 
wall jets, growing on a smooth bed,are presented in Chapter V. 

Experimental findings on the diffusion of a circular wall 
jet in a very wide channel, but growing on a rough bed, are reported 
in Chapter VI. One type of bed roughness was tested and gives useful 
hints on the decay of the maximum velocity along the centerline of 
the jet. 

In Chapter VII, a preliminary investigation on a circular 
wall jet with swirl is presented. The velocity scale is found to 
decay faster when subjected to a swirling motion. Results on the 
growth rates of the swirling jets and the bed shear distribution are 
also presented in the same chapter. 

Based on the present studies, general conclusions are pre- 
sented in Chapter VIII. Some recommendations are then made for 
further studies in this field. 

During the course of this investigation, it was found that a 
similarity analysis can successfully be applied to more complicated 
flows. In the developing zone of a co-flowing circular jet, it leads 


to some simple and interesting results. The analytical results were 
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then verified by experiments, and are presented in Appendix A. 
The present experimental investigations were carried out at 
the hydraulics laboratory of the University of Alberta, Edmonton, 


Canada, during the years 1970 and 1971. 
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CHAPTER I1 
THEORETICAL CONSIDERATIONS 


2.) aelficroduetion 

This chapter presents a theoretical analysis for a particular 
class of three dimensional jets in which the dimensions of the nozzle 
producing the jet are roughly the same in the major and minor axis 
directions. A circular free jet has axial symmetry, but when the 
circular jet diffuses close to a solid boundary the axial symmetry 
is entirely lost, and the flow is three dimensional in nature. As 
the present investigation is primarily on wall jets, any jet of 
circular, square, triangular, elliptic or some other geometrical 
Shape having aspect ratio not very different from unity is referred 
to as a three dimensional wall jet. First a simple rectangular free 
jet is considered for analysis and based on observation of other 
investigators, the Reynolds equations of motion are simplified. These 
simplified or the boundary layer type of equations are then analysed 
to establish the growth rate and velocity decay laws. This approach 
is then extended to the case of three dimensional wall jets. To 
supplement the analysis, and to get the distribution of wall shear, 
a dimensional consideration is presented. As the geometry of the 
flow field does not suggest readily any suitable co-ordinate system, 


the cartesian system is adopted for the sake of simplicity and it 
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also helps to maintain a clear picture of the physical nature of the 


problem. 


2.2 Equations of Motion 
The Reynolds equations in the cartesian coordinates (Schlich- 


ting, 1968) can be written as 
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where x axis is in the axial direction of flow (Figure 2.1), y and z 
axes are normal to the x axis and are along the minor and major axes 


of the nozzle respectively, 


(€.5) 


(8.3) 
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1] 
u,V,wW = the turbulent mean velocity components in the x,y 
and z direction respectively 
u',v',w' = the fluctuating component of the velocities in 
X, y and z directions 


p = mean pressure at the point of consideration 


1S) 
iT] 


mass density of the fluid, and 


kinematic viscosity of the fluid. 


< 
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2.3 Boundary Layer Simplification of the Equations of Motion 


A characteristic feature of turbulent jets is the limited 
extent of the thickness in the transverse direction where the solution 
is sought. Because of the limited extent in the transverse direction, 
jets expanding into infinite regions filled with a quiescent fluid 
nave negligibly small pressure gradients in the x direction. This 
otherwise means that, the constant pressure in the surrounding fluid 
impresses itself on the jet. Another typical feature of submerged 
jet is the smallness of the transverse velocity components, which are 
usually neglected for engineering studies on jets. Based on these 
observations, it is possible to simplify the equations of motion by 
assuming the following. 

(i) Beyond a certain initial reach from the nozzle outlet, 
ye<< xMand ZP<<ext 
(ii) In a major portion of the flow, u >> v and u >>w. Also 
the fluctuating components are much smaller compared to 


the mean values. 
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(iii) The gradients in the y and z directions are much larger 
than the corresponding gradients in the x direction. 
(iv) The pressure gradients in the x, y and z directions are 


all small. 


The desired simplification of tne equations of motion was 
achieved by estimating the order of magnitude of each term in the 
equations of motion (Schlichting, 1968). For this, Eqs. 2.1 to 2.4 
were considered in the dimensional form and the order of magnitudes were 


assigned as, 
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where A is a small quantity, which is at least an order less than unity. 

Since du/ox was of the order of 1, it can be shown from Eq. 
2.4 that dv/oy ~ 1 and ow/dz ~ 1. Furthermore, from the geometry of 
flow, aS y and z become zero on the axis, v and w also become zero. 
Hence in the jet boundary, v ~ A and w ~ A, which was postulated in 
assumption (ii). 

Writing down the order of magnitude of individual terms in 
Eder cob tanougn, Ed. Cac, sue cal DesShnOWethat the entire equation ior 
motion in the y and z directions are of the order of A, and can be 
dropped altogether as an approximation. The turbulent shear stresses 
are usually very large compared to the viscous shear stress and hence 


the viscous terms were neglected. Deleting terms of smaller order 
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in Eq. 2.1, the simplified equations of motion are, 


ace aban pbbe Re, Ue tI aoa fail 
Usacet ay Tins 5y U'vd toxes UW (2.5) 
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These would be referred to as the boundary layer type of equations 

for a three dimensional jet. During this course of analysis, the 

assumption of negligible pressure gradients leads us to the result 
es = 


that, the effect of all turbulent stress components other than u'v 


and u'w' are negligible. 


2.4 Similarity Analysis 


The present studies and those of Trentacoste and Sforza (1967) 
showed that the velocity profiles of the three dimensional jet con- 
tinually broaden. A plot of velocity profile in the physical coordinates 
(presented in Chapter IV) for different x stations and z = 0 plane for 
example, look so very different. This is because, the velocity on the 
axis decreases whereas the profiles get wider. However, interesting 
results are obtained by plotting the ratio of the local velocity to 
the velocity on the jet axis and the ratio of the local distance to 
a characteristic length of the jet. One obvious characteristic length 
of the jet would be the thickness of the jet at the section under con- 
sideration. Exact demarcation of the jet boundary would be a difficult 


task, and alternative length parameters need to be considered. One 


7 i le, 6 ; dest soe 
(2.8) FF Soren 2 shyt v +ihy D tig 
«a . Dae ee bes 1 = 


enattéups to aqyt teyel wisnuuod aid 26 of beristsy sd bivow s2eqT 


ont ,2teyfsns to o21007 ats patind tab (snotensmib sennid aot 


stuesy ond ot av>ebset ana tbsve s1w22estg. sfdtptipan to ‘nakiqwuees 
TVD apit varito ednsnognoa egarte instudys (1s to Jaatts: ait , deny 


-afdipitpan sis "We bre 


ateytonh wivetimte %.5 

(WAC) ssvot2 bns stzonsinsiT to seons bns saneune dnses1q SAT 
-no9 43t. Tenotenanth Seyit aft 70 2alitorg wivotey sit tent bawode 
eeteniiyoos fsotaydq alt at slitorg yisofov to sofq A .nebsord yf fsuntt 
10% sasfq 0 =-s bns anotiste x Fnswsttib 10% {VI +93 qed9 nt batnezsiq) 
sid no vtfoolay srt ,sevsosd ef era? -dnsisttTb yey 02 Aoo! .afamexs 
ontgeovatnt .yevawoH .1sbiw dsp eolftorg sit essionw 2se6ai5eb etxs 
od ytioolsy fso0f oft to oftsy add paittolg yd bentsido o16 ettuest 


od SonstzTb [so00f Shy to: often ait bis etxs Jat. sit no \ivet ev ot . 


dtens! site lyssostans ayotvdo’ snQ.tet afd to iiagnat aid2hvatosyets 6 
-non Yebnw notjoe2 srt 3s 33f ont to eeandotia ‘any od bfuow dst St to 
ffuatVtb 6 od bruow yrsbruad t8t, of Yo notdsovened Joexd .notyersbte 
30 .teisbtano> ed oF bear ape: ntpner shale ee 


14 


such parameter, widely accepted in the study of jets, is the distance 
between the axis and the point at which the velocity is equal to half 
the axial velocity. This is usually referred to as the half width. 
When the velocity profiles at various sections are non- 
dimensionalised by a envoerince choice of velocity and length scales, 
all the profiles collapse into a single curve. This behaviour is what 
we term as "Similarity of the velocity profiles". Thus, in a major 
part of the flow field, the non-dimensional velocity profile is com- 
pletely defined by the non-dimensional distance and is independent 
of the streamwise coordinate x. This is a very important and most 
useful property that is encountered for most of the jet flows. 
For the z = 0 plane, based on the above observations, let 


us assume 


= f,(~) = f (2.7) 
a 'b) (ny) 


where Un is the turbulent mean velocity in the x direction at a distance 
y from the axis of the jet (see Figure 2.1), un is the maximum value 
of ua? bY is a length scale equal to y where Ue Uno! °° ny is non- 
dimensional ordinate and f, represents a functional relation. 

Further in the z - x planes, based again on experimental ob- 
servations (of Sforza and et al. and those to be presented in this 


thesis), we can write 
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where u is the turbulent mean velocity in x direction at a distance 
z from the axis, b, is the length scale in the z or transverse di- 
rection being equal to z where u = u fe (Figure 2.1), ib is non- 
dimensional distance and 9 is a functional relation. Combining 
Eq. 2.7 and 2.8, the local mean velocity at any point can be re- 


written in the following form. 


) + g,(n,) . (2.9) 


Experience on plane and circular jets show that v and w velocity 
components can be non-dimensionalised by the u velocity scale (Anwar 
and Weller, 1969) and the cumbersome operation of getting the dimen- 
sionless v and w components from the continuity equation can be 


avoided. Hence it is assumed that 
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Combining Eq. 2.10 and 2.11 with Eq. 2.7, we can write 
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(2513) 


where ny and Me have been dropped in Eq. 2.12 and 2.13, with the under- 
Standing that f is a function of ny and g is a function of no: 
Similarly in any x-z plane, from experience on self-preserving jets 
(Wygnanski and Fiedler, 1969), the turbulent shear stress components 


can be written as 


y= g4(n,) (2.14) 


and esa (2.15) 
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Combination of these two equations with Eq. 2.7, would give us 
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where p, q) and qo are the unknown exponents. To evaluate these 
exponents, one relationship is obtained from the boundary layer Eq. 
2.5, using the set of assumptions made so far, and following the 


simplifications enumerated below, 
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and the prime on the f function denotes differentiation with respect 


to ny and that on the g functions denotes differentiation with respect 
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From the definition of the similarity of velocity profiles, the left 
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19 
hand side of Eq. 2.29 is independent of x. In order that the equation 
be valid at any section and is nontrivial, the right hand side of Eq. 


2.29 should also be independent of x. So we have, 
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To evaluate the other unknown exponent, we would use the 


integral momentum equation. Since the pressure in the jet is virtually 


. 
OL glgP-f-ght Pd, | Hb. : 
7 


(SE.S) | [=P 70 


20 


constant and equal to the pressure in the surrounding medium, a jet 
issuing into an infinite medium will have the total momentum in the x 
direction preserved and would be independent of the distance from the 


orifice. That is 


4 [ f 6 ucdy dz = 0 (2.34) 
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In Eq. 2.35, the integrand is a finite quantity and is a constant. 


Hence 


od) (2.36) 


That is 


or 2p + q, + 4p = 0. (2.37) 


et 


x 5 boa pane ln “ny 
| | 7; at a0 iebigt - 


a 5 


(8.8) 0 = sb wee || 3 | w 
a, i, ae 
0 = .nb snb dd Se Tt aa 8 1 ib ne 


= 7 

Ai? rey ~ 7 
ee 

- Om 

.inetenoo 6 2t bas -yiktsneup esintt 6 ef basypstat sdf 26.8 pd at 
sons 
| 7 


(a€.8) | ~0= (4 yf 5 q] 4 


(se.8) | 0 = cpt pp + 9S 


~~ ip heh i bis, sien 
a mR NE) 


a) as 
a: ei) 


21 


PVOMeEGSHetc sol and '2asiieap 2a- ell, dsorthat 


b ax (2.38) 


croeevelocity Distribution 


The distribution of velocity for a plane source and axially 
symmetric jet are given by the analytical expressions of Tollmien and 
Goertler (Abramovich, 1963). Tollmien's solutions are based on 
Prandtl's mixing length theory, whereas, Goertler's solutions are 
based on Prandtl's constant eddy viscosity model. The two solutions 
give the velocity distribution with approximately the same accuracy. 
Discrepancies are however observed between the two, near the boundary 
of the jet. 

The analytical solution for three dimensional turbulent 
jets is not well established and the only attempt published in literature 
is that of Sforza (1969). Experimentally it was observed that in the 
far field, the contours describing constant values of streamwise mean 
velocity components in the y-z plane were satisfactorily described 
by ellipses. Also the distribution of streamwise components of the 
mean velocity were similar and identical in various planes, when 


plotted in the nondimensional form. This led to the adaptation 
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of cylindrical coordinates as an appropriate system, and the elliptic 
isovels were transformed to circular isovels by taking the ratio of 

the radius vector r of any isovel to the radius vector ro of a reference 
isovel. The radius vector ro would then be a function of both x and 

the orientation 90. 

The simplified equations of motion were then solved by 
assuming a suitable form for eddy viscosity in terms of ro: General ly 
Speaking, the form of the equations were similar to the one found for 
axisymmetric jets and the approximate solution worked out to be identical 
witn the one given by Schlichting (1968) for the axisymmetric case. 
However, this analysis incorporates the three dimensional jet data in 
determining the constants of integration and the distribution follows 


very closely Goertler's curve (Abramovich, 1963) for a free circular 


jet. 


2.6 Extension to Three Dimensional Wall Jets 
Gr Orme rec Us Olav 1SCOS LY 
In the presence of a solid boundary, the molecular viscosity 
cannot be neglected completely. Near the wall the velocity gradients 
are large, and so, for a low-viscosity fluid flowing past a solid 
boundary, the flow field can be divided into two regions: a comparatively 
thin layer of fluid located close to the solid wall where viscosity is 
important and the outer part in which viscous effects can be neglected. 
Glauert (1956) considered the laminar and turbulent, radial 


and plane wall jets, and showed that a complete similarity in velocity 
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profile in the wall region and the outer free turbulence region is 
not feasible as a result of the viscosity model being different in the 
two regions. It was also shown that the concept of constancy of the 
flux of exterior momentum flux is not applicable for a turbulent wall] 
jet. Glauert's simple power laws, which express the local rates of 
change with x of maximum velocity and the jet width, have the exponents 
themselves varying slowly with x. In later experiments (Schwarz and 
Cosart, 1960), it was found that a single velocity scale and a single 
length scale correlate all the velocity data; and leads to the concept 
that it is possible to have experimental similarity, if not theoretical 
similarity, over tne entire turbulent thickness of the jet. The 
equation of motion indicate the single length scale to vary as yO 
and the velocity scale to vary as xP The exponent p was determined 
empirically, since a second relationship could not be found out for 
the plane turbulent wall jet. 

In view of the difficulties encountered, to account for the 
viscosity in a plane wall jet, no simple approach is at sight for a 
three dimensional case. In the hope of establishing the exponents 
which might serve as a good approximation and provide reasonable 
answers for engineering calculations, the shear stress on the bed of 
the channel is neglected. This expectation seems to be borne out in 


the light of agreement between the predicted laws and the experimental 


results. 
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2.6.2 Simplification of the Equations of Motion 
Following the line of argument presented for a free jet, 
the equations of motion (Eq. 2.1 through Eq. 2.4) for a three dimensional 


wall jet can be simplified as 
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2.6.3 Similarity Analysis 

The non-dimensional velocity distributions are represented 
by the forms given in Eq. 2.7 and Eq. 2.8, with an exception regarding 
the definition of the length scale b,- For wall jet studies the length 
scale bY equals y, the height measured from the wall, where Ti aig Uno! < 
and du, / oy is negative (Figure 2.2). All other assumptions and forms 
of representation remain identical with the free jet case. 


Subs.ti-tutingsforathe terms in Eq. 2.39, 1t can be shown, as 


in Section 2.4, that the following expressions will be obtained. 
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For turbulent jets UnoPy!¥ is generally large, so that the first 

term on the right side could be neglected and this assumption would 
run into trouble only close to the wall in the laminar sublayer. 
Further for a wall jet b/d, is less than unity, so that the second 
term on the right hand side of Eq. 2.41 can as well be dropped. For 
Similarity, since the left hand side of Eq. 2.41 is independent of x, 
the right hand side should also be independent of x. This would give 
the exponents the Gls © 1, the details being worked out in Section 


(baht es 


2.6.4 Momentum Integral Equation 


First of all, the momentum integral equation is worked out 
taking into account the wall shear. The loss in momentum because 
of the wall shear is estimated, which then leads to a simplified 
momentum integral equation that successfully gives the exponent p. 


From Eqs. 2.39 and 2.40, we can write 
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Equations 2.42 is integrated over the space bounded by y = 0 plane 


and reduces to 
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where To is the shear stress on the boundary. To evaluate the integrand 


on the right hand side of Eq. 2.43, it is rewritten as 
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x' = streamwise distance measured from the nozzle section 
Perini boundary shear stress on the centerline and is in 
the x direction. 
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and r | 

we = tay (2.47) 
0 (x'/VA) 
Here, ays Ce and C are empirical constants, provided from experiment. 
U5 is the efflux velocity at the nozzle outlet and A is the area of 
cross section of the outlet. The form of Eq. 2.47 is a simplified 
version that fits the experimental data with no high degree of accuracy. 
More elaborate expressions are presented for the maximum velocity 
decay in Eq. 4.1. However, the final answers for momentum loss do 


not differ very much even when the more precise expressions are used. 


Making use of Eqs. 2.45 to 2.47 in Eq. 2.44, we have 
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where, Mo = pAUy” = efflux momentum of the nozzle. 
Integration of Eq. 2.48, would give 
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The constant of integration a5 is found from the condition 
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that at the virtual origin x' = Xo» and the integrand on the left 
hand side of Eq. 2.49 is equal to the efflux momentum flux Mo: Hence, 


Eq. 2.49 can be written as 
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Based on the present experimental data (discussed in Chapter IV), the 


empirical constants involved in the above expression are given as 


C = 8.4 (approximate value) 


0.206 


iT} 


cy 


ay 0.0065 


Equation 2.50 is then simplified to 


ik | 0 uc dy dz = My (1 - 0.0476 Xn a) : W250) 
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Noting that XQ iS approximately ten times the nozzle thickness, the 
loss in momentum in a distance of one hundred times the nozzle thick- 
ness is around 11% of the efflux value. This is very close to the 
value estimated for the case of a plane turbulent wall jet (Newman, 
1969). This demonstrates that as a rough approximation the three 
dimensional wall jet momentum is conserved and the integral equation 


(Eq. 2.51) simplifies to 
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From Eq. 2.52 and the values of qy and qo determined in 
Section 2.6.3, the exponent p = - 1. 


hence, for a three dimensional wall jet, we have 


be ox (22.53) 


For three dimensional wall jets no analytical solution for the distri- 
bution of velocity is available at the present time and is established 


from experiments. 


<./ Dimensional Consideration 

The velocity decay and the growth of the jet can conveniently 
be predicted from simple dimensional considerations. Experience with 
the plane jet (Newman, 1961) diffusing in an stationary ambient show 
that the efflux velocity and the thickness of the slot do not affect 
the problem individually. The two independent variables can be re- 
placed by one variable, the efflux momentum flux My: This leads to 
the hypothesis that even for a three dimensional jet, the one variable 
Mo would be a dominant parameter. 

In most practical cases the Reynolds number is greater than 


4 


10° and it would be reasonable to neglect the effect of viscosity. 


Hence, the maximum mean velocity along the centerline would be determined 
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by the momentum flux at the nozzle section, the fluid density and the 
streamwise distance x measured from the virtual origin (discussed in 


detail in Chapter IV). 


i.e. Usa (Mp 50 sx) ; (2.54) 


Using the Pi theorem, Eq. 2.54 can be written as 


HO age (2.55) 
/My/ox® 
where, C iS an empirical constant. 
Expressing the momentum flux in terms of A, the area of cross-section 
of the nozzle, and the efflux velocity Up» Eq. 2.55 may be rewritten 


in the following form 


U 
1k See (2.56) 


The form of Eq. 2.56 looks very general, and suggests that this would be 
valid for any arbitrary shape as long as the aspect ratio of the 
nozzle does not differ very greatly from unity. 


For a circular jet of diameter d, Eq. 2.56 becomes 
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and for a square orifice of side B, 


“mo fe 


where, C, and C, are empirical constants. 
On similar lines of argument, the problem regarding the 


length scales can be stated as 


by» by = $(My.0 x) (2.59) 

and would lead to 
by = C3x (2.60) 
b, = C 4x (2.61) 


Here again, C, and Cy are empirical constants. For a three dimensional 
wall jet, the effect of the viscosity is predominant in the laminar sub- 
layer close to the wall and may not affect the relationship 2.56, 2.60 
and 2.61 to any great extent. However, the shear stress on the wall 

is a result of the viscosity of the fluid, and should therefore be 
considered. The shear stress in the center-plane can then be expressed 


as 


Tom = (Mg 2? »Xsv) (202) 
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shear stress on the wall for z = 0 values, 


Ss 
iT} 


kinematic viscosity of the fluid. 


Using dimensional considerations, Eq. 2.62 can be reduced 


to a general form 


T 2 UWA 
om (y= F(A (2.63) 
pUy VA 
Further, Eq. 2.63 may be represented as 
G 2 U.d 
re (q) es cae) (2.64) 
pu, 
for the circular nozzle, and 
af c UB 
om (A) = (SY (2.65) 
PU 


for the square nozzle. 


Summary 


An order of magnitude study is presented to simplify the 
equations of motion. Based on the observations that the velocity 
profiles in the x-y and x-z planes are similar, the velocity decay and 
growth rate laws for a three dimensional free jet are established in 
simple exponential forms. From experimentally derived constants, an 


estimate for the loss in momentum flux of a three dimensional wall jet 
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is presented showing that there would be a loss of about 11% in a distance 
of 100 times the nozzle size. Therefore, as an approximation, the 
momentum of a three dimensional wall jet is assumed to be preserved. 

Based on this approximation, the velocity decay and growth rate laws 


are seen to be identical with the free jet case and are given as 


A dimensional consideration for the flow is also presented taking 
into account that the outlet velocity and the area of cross-section 
of the jet do not appear as two separate variables, and rather combine 
to give a single variable-the efflux momentum flux. The variation of 
the velocity scale and the length scales, found from dimensional 
considerations, are in agreement with the relationships established 
from the equations of motion. A form for the variation of the bed 
shear stress along the centerline of the jet is also presented. 

To check the predicted relationships, some of the available 
data on three dimensional free jets are analysed in Appendix B. The 
wall jet investigations are, however, presented in great detail in the 


subsequent chapters. 
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CHAPTER III 
EXPERIMENTAL ARRANGEMENT AND 


MEASUREMENT TECHNIQUES 


Sele introduction 

This chapter describes the experimental arrangements adopted 
and the limitation of the measurements. A set of preliminary experi- 
ments were conducted in a readily available 12.5" wide flume and an 
extensive investigation was later carried out in a wider flume. The 
narrow flume would be designated as flume 1] and the wide one as flume 2. 
Effects of the bed roughness, side walls and the swirl were all studied 
in flume 2. Probes for the measurement of bed shear and the velocity 
vector in three dimensional flow field had to be built and calibrated 
when standard calibration charts were not available. The swirl 
generating mechanism, used in the present experiments, is also dis- 


cussed in detail. 


3.2 Experimental Flume 1 

This flume is 12.5" wide, 18" deep and 15'-6" long, with a 
horizontal aluminum bottom and plexiglass sides. Figure 3.1 gives a 
schematic representation of the flume. There is a flap gate at the 
downstream end of the channel which controls the tailwater depth, and 


discharges into the laboratory sump. A constant head was maintained 
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SECTIONAL VIEW (not to scale) 


FIGURE 3.1 EXPERIMENTAL FLUME 1 (SCHEMATIC REPRESENTATION) 


38UT AIH JAIOT =) at } 


ey Sa an a ay ed we ee aes oa teen | om Ss i 
: 7 a re 


| i 
HAT 
WAS 


LAV Ve SSS SSTASSSS SSS SS SS = S35 we 


oo 
i 

s : 
: 
aA 
4 on) Er 
{elese of ton) W3IV JAAOITIZ2 : 
i ’ 

) 

i _ OW 


! (MOUTATMERARIEA OLTAMAND2) J SMULY JATMOMIARIKA £.e 


37 
at the nozzle by a vertical stand-pipe 6" diameter and about 11'-6" high 
having an overflow device. The stand-pipe was fed by the city water 
supply. Near the bottom, the stand-pipe is turned through 90° by an 
elbow, and connected to a contraction. The contraction was prepared 
by joining two commercially available reducers of 6" x 2" and 2" x 1" 
Size and the flow was contracted in a short length of 8". Test nozzles 
of the desired shapes were turned from solid brass cylinders 1" in 
diameter and of a short length of 2" to keep the boundary layer growth 
to a minimum. As far as possible, bell-mouth entries were provided in 
the nozzle and were fixed in the contraction without any leak. The 
test nozzle that fits the contraction, rests symmetrically at the up- 
stream end of the channel and after some trials the nozzle was located 


tangentially to the bed. 


3.3 Experimental Flume 2 
This is a rectangular tank 43" wide, 48" deep and 11'-4" 


long with plexiglass sides. To conduct wall jet experiments, a false 
bed was provided at a depth of 26.5" from the top. The smooth aluminum 
plate forming this bed was 43" wide and 9'-0" long, resting on rigid 
channel sections and a system of jacks that allowed to level the bed 

as closely as possible (Figure 3.2). The flow in the system was re- 
circulated by a centrifugal pump having its intake in the downstream 
end wall of the tank. A constant head at the test nozzle was main- 
tained by an arrangement similar to the one described in flume 1. 


However, to maintain a constant depth of flow in the flume during 
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the tests, the overflow pipe discharges back into the experimental 
tank behind a baffle wall. The contraction arrangement and the details 


of the nozzle were in general identical with those in flume 1. 


3.4 Roughness Studied 


A continuous strip of #36 aluminum oxide wet-or-dry cloth 
manufactured by the 3M company was glued to the smooth boundary. The 
grains of this cloth were very angular and are shown in Figure 3.3. 
The grains were imbedded in a water-proof glue to varying depths. 

The grain size ranges from 0.38 mm to 0.70 mm, with a median size of 
0.56 mm(0.00183 ft.). This was taken as the nominal roughness height 
K. 

While measuring bed shear on rough boundaries , Hollingshead 
(1972) found that the equivalent sand grain roughness for this parti- 
cular cloth is 3.3 times the nominal height (i.e. K. = 0.0061 fim) 

A high ratio of K./K = 3.3 resulted as the grains were very angular 
and many grains that had their long axis about twice the sieve size 


were oriented with the axis normal to the boundary. 


3.5 Velocity Measurements 
3.5.1 Total Head Tube 


The velocity field for jets on smooth wall and no swirl was 
explored by a total head tube 0.05" in external diameter and 0.033" 
internal diameter. During some of the tests in flume 1, the static 
pressures were measured by a Pitot-static probe and the distribution 


was essentially found to be hydrostatic. Therefore, in most of the 
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FIGURE 3.3 SAND GRAIN ROUGHNESS 
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4] 
experiments a static pressure tapping was located on the side wall of 
the flume and was assumed to indicate the static pressure throughout 
the jet thickness. For lower range of velocities, the total head tube 
and the static pressure tapping were connected to a manometer board 
that was inclined at a slope of 1:2 and read to the nearest + 0.05" 
of water. For higher range of velocities, the two pressure leads 
were connected to a U tube mercury manometer and the pressure dif- 
ferential was measured with an accuracy of + 0.1" of mercury. 

During the investigations on the effect of side walls, the 
water surface was not level, and the static pressure had to be 
sensed by a Pitot-static probe stationed close to the total head tube. 
In all cases, the probe was mounted on a traversing mechanism that 
could be located to the nearest 1/1000 of a foot in the transverse 
and vertical direction, and to the nearest 1/8 of an inch in the 
streamwise direction. The corrections to account for the turbulent 
fluctuations, tne viscous effects and the displacement of the effective 
center were generally believed to be small and no correction was 
applied to the measured velocity. The reading of a Pitot tube is 
very insensitive to yaw and an angular displacement of 20° would 
affect the reading by less than 0.5 percent of 0 u=/2 (Rosenhead, 1963). 
Hence, in a major part of the jet, the correction for angle of attack 
would be very small and no attempt is made to correct the present 


velocity measurements. 
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3.5.2 Fixed Direction Five Hole Probe 

For velocity measurements in swirling jets, the conventional 
Pitot tube would not suffice as the angles of attack will be very 
large. In addition the present investigation was aimed at measuring 
all the three components of a velocity vector which would require the 
measurement of the yaw and pitch angles. In order to find the unknown 
quantities - the static pressure, the total pressure and the two angles 
of inclination of the velocity vector with a reference direction, a 
five hole probe was used. The five hole probe can be used as a null 
probe by rotating it, till its principal plane coincides with one of 
the planes containing the velocity vector. This type of operation 
was not possible in all cases, and insensitivity of the probe to 
angles of attack combined with the slow response of the manometers, 
particularly for small size tubes, makes it hard to align the probe 
sccuracely with the velocity vector. Hence,a fixed direction five 
hole probe was preferred inspite of its complex calibration curves. 

In the present case the probe was prepared from stainless 
steel tubes 0.05" outside diameter and 0.033" inside diameter and 
has a 45° chamfer. Lewis (1965) has given calibration curves for a 
hemispherical nosed five hole probe, but the curves are not standardised 
and any change in the nose shape would change the calibration curves 
significantly. Hence, the present probe was calibrated and the principle 
involved is outlined below, following in general the line of thought 
adopted by Lewis. Let Py> Po» etc. be the pressures recorded by the 
five holes (numbered as shown in Figure 3.4) and Po be the static 


pressure. Then 
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Be = Py Ky oa i = 1,2,3,4,5 (eal) 


where q is the dynamic pressure. 


AXIS OF PITCH HOLES 


STAINLESS STEEL TUBE 
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FIGURE 3.4 DEFINITION SKETCH FOR PITCH AND YAW 


Using Eq. 3.1, the following coefficients were defined 
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It was assumed that all these coefficients were only func- 
tions of the pitch angle ¢ and the yaw angle w, and can be represented 
in an unique manner. 

The probe was calibrated in the potential core of a plane 
wall jet having a velocity of 6 ft./sec. The dynamic pressure was 
measured by a standard Pitot-static probe. During calibration, values 
of w up to + 55° could be attained, but the test facility restricted 
the values of ¢ to a maximum of 45°. Both o and w could be set to 
the nearest degree and were varied at intervals of 5°. For each com- 
bination of » and wy, the pressure sensed by all the five holes were 
recorded. Because of the restrictions on setting negative > angles, 


it was assumed that the probe is symmetrical about a plane passing 
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through the centers of tube 2, 1 and 4. Thus, in the typical calibra- 
tion charts presented in Figure 3.5, Ke for negative » iS a mirror 
image of the values for positive ¢ about the zero pitch axis. How- 
ever, for the Ka coefficient the values for negative 9 were the reci- 
procals of the corresponding quantities for positive » angles. The 
other two coefficients Kyo and K were symmetrical about the o = 0 
axis (Figures 3.6 and 3.7). 

At the point of measurement the five pressure readings are 
known, and K, and Ke are computed. From the calibration curve 
(Figure 3.5) the corresponding » and ~ angles were noted. From 
Figure 3.6, Ky2 value for the noted angles is read out and the dynamic 
pressure q found from Eq. 3.2. The value of q should work out to be 
a positive quantity. A negative q would indicate that there are other 
possible pairs of wy and ¢$, that correspond to the known K, and Kes 
needing a new estimate. K, is then read from Figure 3.7 and with Py 
and q known, Pp can easily be computed. 

In certain parts of the chart (Figure 3.5) the coefficient 
K 


and K, become infinity and makes it impossible to estimate the 


7 8 
angles ¢ and y. However, this can be overcome by plotting charts 

for the other coefficients and using the one that gives finite values 
for reliable interpolation. 


The three velocity components in terms of ¢ and y are given 


by 
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FIGURE 3.5 CALIBRATION 
FACTORS K7 AND Kg FOR FIVE HOLE PROBE 
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u = V(1+tan@u+tan 


oa ee 
v = V tano(1t+tanW+tan* ) 3174) 


Lip ieee Mle 
w = V tanv(1+tan-yt+tan“9¢) 


The probe was calibrated in an uniform velocity field away from the 
wall, whereas the conditions during measurement happen to be quite 
different. The probe was used in a velocity field having large 
gradients and would introduce errors. This was particularly acute 
in a field close to the nozzle section. The viscous effects and the 
proximity of the wall would introduce additional errors in measure- 
ment and there seems to be no information available to account for 
this. There is also a limitation on the smallest size of probe that 
can be made from the available tubes. In view of the complex nature 
of the calibration curves, no comment was made on the order of errors 
involved. It is believed that the probe can measure @ and y to an 
accuracy of + 2° and the u component can be measured very reliably. 
The accuracy in the measurement of small v and w components seems 


to be limited. 


3.6 Boundary Shear Stress Measurements 
3.6.1 Preston Tube 


The Preston tube is a total head tube resting on a boundary 
and held against the flow. This is a convenient method of measuring 


the boundary shear when the direction of the flow is known. The 
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stagnation pressure was found from the sensing device connected to the 
total head tube, and the static pressure was found by different means. 
The difference between the two pressures Ap is correlated with the 
boundary shear stress Tg: The correlation was developed by Preston 
(1954), mainly from similarity considerations and the existance of 


the law of the wall. This correlation is given in the form 


nN q@ inde 
=e = F(—s) (3.5) 
2 2 
4ov 4ov 
where Ap = dynamic pressure 

d = outer diameter of the total head tube 

Oo = mass density of the fluid 

vy = kinematic viscosity of the fluid 


The principle of the method has been confirmed by Head and Rechenberg 
(1962), and Rajaratnam (1965) derived calibration curves theoretically 
by assuming a suitable velocity distribution in the wall region. The 
Preston tube technique is found to be valid even in regions with 
milder pressure gradients and have been investigated extensively by 
Patel (1965). In the present investigation, Patel's curves have 

been adopted. The calibration curves available are for cylindrical 
Preston tubes with square ends, having a ratio of the internal to the 
external diameter of about 0.6. Even though this ratio is not very 
critical, to use the available curves with high accuracy, the diameter 


ratio should be between 0.5 and 0.7, preferably equal to 0.6. 
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In the present measurements, the total head tube used for 
taking velocity profiles was also serving as the Preston tube and had 
0.05" outside diameter and 0.033" inside diameter. The static pressure 
tap was located on the side wall of the flume. Temperature of water 


9 ft.*/sec. 


was around 66°F and v = 1.11 x 10° 
The method requires that the size of the tube be small com- 
pared to the boundary layer thickness so that it would lie in a zone 
where the law of the wall is valid. However, in many of the present 
measurements the tube was larger than the wall law region and in some 
instances the tube occupied around 35% of the boundary layer thickness, 
casting some doubt on the accuracy of the boundary shear results. 
Sigalla's (1958) investigations on the effect of instrument size on 
the measurement of skin friction in a boundary layer showed that the 
results would be within 4% accuracy, even when the diameter of the 


tube was about 40% of the boundary layer thickness, and therefore, the 


present measurements were acceptable. 


3.6.2 Yaw Probe as a Preston Tube 

When the direction of shear stress is at some unknown angle 
with the free-stream, in principle its magnitude and direction can 
still be obtained by rotating the Preston tube till the maximum dynamic 
pressure is indicated. Since the dynamic pressure recorded was very 
insensitive to the angle of yaw 0, particularly when the angle was 
less than 20°, this approach is very limited. Further, installing a 
mechanism to rotate the tube would be an added inconvenience and time 


consumming. In order to overcome these difficulties, Rajaratnam and 
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Muralidhar (1967) used a yaw probe on the boundary and gave the necessary 
calibration curves for shear computation. The nose of the probes was 
chamfered at 45°. The principle involved in this method of shear 
measurement is described briefly in the following paragraphs. 

The yaw probe was placed on the boundary facing in the 

reference direction, the x axis. If the shear stress was at an angle 
8 with the reference direction and if the center tube was designated 
as number 2, with the tube on the side of shear stress as number 3 
and the remaining tube as number 1, the total pressure recorded by 


them would be, 


Ap (3.6) 


where Kig> Kage K39 are the respective calibration factors and are 
only functions of 8, Po is the static pressure and Ap is the dynamic 
pressure. Combining the various expressions in Eq. 3.6 and 3.5 the 


following can be written 
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By suitably combining the expressions in Eq. 3.7, one can 
write 


Ka,-K 


UB 2pepec0a COP ag, 


2 | (3.8) 
Pi-Po = Kyg7kog «=O 


The above expression indicates that Ky is only a function of 9 and 


can be conviently calibrated. (Figure 3.8). 
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FIGURE 3.8 CALIBRATION FACTOR Ky FOR YAW PROBE 


(Rajaratnam and Muralidhar, 1968) 
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From Eq. 3.6, it is possible to write, 


(379) 


In measuring shear stress at any point, Py» Po and P. are known 
quantities so that Ky can be computed. From Figure 3.8 the corres- 


ponding 8, and then Ki and K59 from Figure 3.9 were noted. 
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FIGURE 3.9 CALIBRATION FACTORS FOR YAW PROBE ON THE BOUNDARY 


(Rajaratnam and Muralidhar, 1968) 
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59 
Since all the quantities on the left hand side of Eq. 3.9 were known, 
ae was determined. Using Patel's calibration curves ce and hence 
Tg were computed. 

It is to be noted that the law of the wall for a three di- 
mensional boundary layer seems to be valid in a much smaller region 
than the corresponding two dimensional case, and therefore, the tubes 
should be as small as possible. Stainless steel tubes having an out- 
side diameter of 0.045" and an inside diameter of about 0.6 times the 
outside were used in the present investigation. The nose shape was 
made geometrically similar to the one calibrated by Rajaratnam and 


Muralidhar (1967) and hence their curves have been used. 


3.6.3 Preston Tube Technique Applied to Rough Walls 


Hwang and Laursen (1963) extended the Preston tube technique 
to measure shear on rough boundaries. A relationship between the dynamic 
pressure acting on a Pitot tube in contact with a rough boundary and 
the local boundary shear was worked out analytically using logarithmic 
velocity distribution with Nikuradse's constants. For the fully developed 


rough flows the final relationship worked out as, 


2 2 4 
Ap 2°46.531 {Llog wo]. - 10g 210.25 (2) + 0.0833 (2) + 1] 
T K K h h 
0 S S 
a 2 a 4 
+ [0.25 (7) + 0.1146 (7) + |} chery) 
where h = height of center of stagnation tube from zero datum 


inner radius of the stagnation tube 


e3) 
iT} 


K. = equivalent sand roughness. 
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The above equation converges rapidly for small values of a/h. The 
pressure to shear ratio which depends on the roughness size a/K. and 


the relative position of the tube (h-a)/K. is indicated in Figure 3.10. 
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FIGURE 3.10 PRESSURE SHEAR RATIO FOR ROUGH BOUNDARIES (ANALYTICAL) 


(Hwang and Laursen, 1963) 
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For measurement in the hydraulically smooth and transition 
region (u,K/v < 70) empirical curves for applying correction are given. 

In the present measurements a stagnation tube having outside 
diameter 0.11" and inside diameter 0.085" was used. The boundary had 
an equivalent roughness eo 0.0061 ft. Assuming Z9/K, = 0.15, where 


z, is the height of bottom of stagnation tube from zero datum, 


0 
(h-a)/K. =9).33 and) d/ K-agUs0/ Jveehhom blgureso. 10, 


ae eye ea) 
To 


This expression was used to compute the boundary shear in hydraulically 


rough flows. 


3.6.4 Velocity Profile Method 


For a two dimensional boundary layer, the wall shear stress 
was assumed to be very nearly equal to the turbulent shear outside the 
laminar sublayer. Further, on the assumption that this shear stress 
remains constant in a small zone close to the wall, an expression for 


the velocity distribution can be derived as (Schlichting, 1968) 
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rough boundaries (i.e. u,k./v > 70), B, = 8.5. Regardless of whether 
the flow is hydraulically smooth, transitional or rough, adopting a 


value of 0.4 for k, the above equation can be rewritten as 


T 
/2= uy, = are (3.13) 


iy 
log e 

This is a convenient expression for boundary shear computation 
as a knowledge of K. and B. are not needed. 

To use this method, it would be necessary to determine the 
correct datum from which y is to be measured. The datum selection was 
generally made by trial and error by plotting u versus log y, y being 
measured from an assumed datum, till a straight line plot was obtained. 

The displacement of the effective center of the Pitot tube 
was another important aspect that needed consideration, particularly 
for measurements taken very close to the wall. The exact correction 
applicable for a Pitot tube resting on a rough boundary is open to 


question, although for a smooth boundary it is around 0.18 times the 


tube diameter (Macmillan, 1957). 


3.7 Swirl Generation 
Se) lee SWI) Genera Cor, 

The experimental layout for producing swirling jets is sche- 
matically represented in Figure 3.11. The generator had an outer 
casing 6" in diameter and an inner steel pipe having a inside diameter 
of 1" and 0.5" thick wall. The inner pipe was 18" long and had four 


tangentially cut slots 0.19" thick extending for a length of 4.5". 
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FIGURE 3.11 SWIRL GENERATOR (SCHEMATIC REPRESENTATION) 


At the exit end of the inner pipe, a 0.75" nozzle having rounded entrance 
was provided. Quantities of water let into the axial and the tangential 
supply pipe could be varied. The highest swirl was produced when the 
entire flow was led through the tangential slots, whereas, lower swirls 
could be produced by diverting only a part of the flow through the 

slots. However, in all cases the stand pipe was made to overflow so as 


to maintain a constant head. 


3.7.2 Discharge Measurement 


The discharge rate through the tangential slots as well as 
the total discharge in the main line were measured by sharp edged 


orifice meters 1" in diameter, with pressure taps located at 1.5" up- 
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stream and 1.0" downstream of the orifice plate. The discharge was 


computed from the flow equation (Spink, 1958) 


eA 
Q, = oND hake Gy v bie (3.14) 
where Q = rate of flow, units same as N 


a = an operating figure, depending on ratio of orifice 
and pipe diameter 

N = a constant which takes care of units of measurement 

D = diameter of the pipe in inches 

F_ = ratio of area of primary device opening at flow 
temperature to that at 60°F 


F_ = manometer factor for liquid on mercury 


m 

Gy = specific gravity of liquid at flowing temperatures 
G, = specific gravity of liquid at base temperature 

hy = operating differential as read on the differential 


gauge in inches of water 


In the present setup D = 1.5", hence a = 0.3074, N = 0.01263, E =i), 
VG -/Go = ] and Eq. 3.14 can be simplified to 


Q, = 0.00875 vAH cfs ... (3.15) 


where AH = differential manometric head in inches of water. 
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3.7.3 Measurement of Swirl Number 
For a circular jet, dimensional analysis (Chigier and 
Chervinsky, 1967) and integration of the equations of motion show that, 


the degree of the swirl can be characterised by a single parameter, 


t 


S = (8816) 
(P+M) rq 
where S = swirl number 
T = angular momentum of the jet 
P = total pressure in the jet 
M = axial momentum of the jet 


aye radius of the nozzle 


The form of Eq. 3.16 is based on the fact that the quantities T and 


(P+M) are preserved for a free jet. Further ijt can be shown that 


T= | 0 réu w dr (Sei7) 
O 
and = 
2 we 
P+M = | o(u aio) Y dh ac C3e1c) 
0 


The swirl number defined in Eq. 3.16 can be varied by changing the 

ratio of the discharge through the tangential slots to the discharge 
in the axial direction of the generator. This relationship was pri- 
marily guided by the geometry of the slot and the size of other com- 


ponents of the chamber. There seems to be no way of establishing the 
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Swirl number that could be generated for different discharge ratios, 

except by direct measurements. The ratio was set at four different 

values in the range of 0 to 1 and for each ratio, velocity profiles 

in a radial plane of the free jet was measured by the five hole probe. 

From the measured values of u, w and r, the quantities on the right 

hand side of Eqs. 3.17 and 3.18 were computed graphically. For a particular 
discharge ratio, the measurements were conducted at distances of 2" 

and 4" from the nozzle section. Theoretically the swirl number should 

be the same for any section, but the present results showed some vari- 


ation (Table III-1) and is perhaps because of experimental errors. 


TABLE III-1 
SWIRL NUMBER MEASUREMENTS 


Discharge ratio Swirl No. S$ Average 
Qa/07 = 3220, Xie AeOE S 
Sweet 
jet) Ui) Uzz200 D200 
0.877 0.205 02197 e201 
Uno 0.10 mae 0.10 
G2a05 Q707c5 070535 0.063 


Q_ = discharge through the slots 


Qi = total discharge 
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The average values of the swirl numbers are plotted in Figure 3.12 as 
a function of the discharge ratio and serves as a calibration curve 


for the swirl generator used in the present investigation. 


FIGURE 3,12 SWIRL NUMBER-DISCHARGE RATIO RELATIONSHIP 
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CHAPTER IV 
THREE DIMENSIONAL WALL JETS 


4.1 Introduction 

Wall jets of circular, rectangular, triangular or other 
Similar geometrical shape, having finite aspect ratio are denoted as 
three dimensional wall jets. Experimental investigations on three di- 
mensional wall jets in wide channels and on a smooth wall (or bed) 
having no swirl is presented in this chapter. The other experimental 
investigations are those of Sforza and his associates at the Polytechnic 
Institute of Brooklyn (1966, 1967) and are limited to rectangular out- 
lets. Based on present observations, relationships for the decay of 
maximum velocity, growth of the jet and bed shear distribution are 
presented in suitable forms. Experimentally found relationships are 


then compared with theoretically established laws. 


4.2 Review of Existing Work 


The bulk properties of turbulent, incompressible, rectangular 
wall jets with air as the media, issuing into a quiescent ambient have 
been experimentally studied by Viets and Sforza (1966) and Sforza and 
Herbst (1967). Depending on how the maximum velocity decays, the flow 
field was divided into three distinct regions - the potential core, 
the characteristic decay region and the radial type decay region. 

(a) Potential core region (P.C.): The maximum velocity remains 


almost constant and is close to the exit velocity. This 
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65 
region extends from the nozzle section to a distance of 
about ten times the nozzle height. 

(b) Characteristic decay region (C.D.): The maximum velocity 
decays as a constant power of x. Although the mixing from 
the near boundaries of the outlet has reached the center of 
the flow region, the mixing from the far boundaries that 
are on the sides have not extended over the entire flow 
field. Hence,in this region the velocity decay was re- 
lated to the aspect ratio of the nozzle. The centerline 
velocity profiles normal to the wall was found to be similar. 

(c) Radial type decay region (R.D.): In this region, the velocity 


es and the mixing from all boundaries of the 


decays aS xX 
nozzle extend over the entire flow field. The flow tends 

to loose memory about the nozzle geometry. The centerline 
velocity profile taken normal to the wall and the profiles 


in the transverse direction were found to be similar. 


Other important features noticed were that the growth of the 
mixing layer normal to the wall was independent of the orifice shape. 
Further, the near-field growth in the transverse direction was found 
to be affected by the nozzle geometry and the transverse growth rate 
happened to be larger than the growth in the normal direction. 

Irregularities in velocity profile of the C.D. region were 
observed in transverse planes and was believed to be caused by the 
induced velocity connected with a system of vortex rings surrounding 


the jet. 
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4.3 Experiments 

Two series of experiments were conducted. The first series 
of four experiments were conducted in flume 1 and the second series 
of six were conducted in flume 2. Descriptions of the two flumes are 
given in Chapter III. In the second series, the fifth run was for a 
free square jet and the results are presented in Appendix B. The 
shapes studied include a circle, a square, an ellipse, an equilateral 
triangle with its base resting on the bed and two rectangles. The 
significant details are given in Table IV-1. 

In all experiments, the nozzle was deeply submerged so that 
the water surface was level thereby ensuring a negligible longitudinal 
pressure gradient on the wall jet. The velocity field was explored 
with a total head probe 0.05" external diameter and in most cases, the 
static pressures were obtained from a side wall tap. The boundary 
shear along the centerline of the jet was explored by the same total 
head tube, used as a Preston tube. The off-center shear stress dis- 


tributions were,however, measured by a yaw probe 0.045" thick. 


4.4 Typical Experimental Results 


A typical distribution of velocity (u_) in the center-plane 
is shown in Figure 4.1, where the distance x' is measured from the 
nozzle section. It is seen that the velocity distribution assumes 
the wall jet shape very quickly. Figure 4.2 shows another typical 
plot of variation of velocity u with z in different horizontal planes 
for a number of x stations. In any one plane, the maximum velocity 


that occurs on the centerline decays as the distance from the nozzle 
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FIGURE 4.1 VELOCITY DISTRIBUTION IN THE CENT ER -PLANE 
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is increased and simultaneously there is a continuous widening of the 
profile. The distribution of wall shear is given in Figure 4.3 and 
its behaviour is in general similar to the transverse velocity profile. 
The boundary shear stress along the centerline is found to decrease 


with x and some of the significant results are tabulated in Appendix C. 


4.5 Analysis of Experimental Results 
4.5.1 Velocity Distribution in the Center-Plane 


For all the shapes tested, the velocity profiles in the center- 
plane were checked for similarity by plotting rie le versus hy: A typical 
plot for the run #8 is shown in Figure 4.4. It was found that for every 
case Studied, the velocity distributions are similar after a streamwise 
distance of about 15 times the nozzle height. The rest of tne non- 
dimensional velocity distribution for other shapes are given in Appendix 
C. The curve of the plane turbulent wall jet on smooth walls with zero 
pressure gradient, denoted for the sake of brevity as the classical 
wall jet (CWJ), is plotted in Figure 4.4 for comparison. Although 
the experimental points could be well compared with Glauert's (1956) 
analytical solution by selecting a proper value of a, it is thought 
that a proper selection of a is somewhat arbitrary and hence the CWd 
curve derived from most available two dimensional wall jet experimental 
data would be more reliable. The classical wall jet curve describes 
the present data reasonably well and there appears to be some scatter 
for large values of ny and is mostly because of the large errors in- 
volved in measuring small velocities. There is also scatter in the 


data for ny less than 0.1, indicating that a single length scale 
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is not adequate. For large x, this scatter gets small and for all 
practical purposes the CWJ curve is satisfactory over the entire 
thickness. Further it is seen that in all cases, the maximum velo- 


city occurs at about 0.16 times by from the wall. 


4.5.2 Velocity Distribution in the Transverse Planes 


Consider the velocity profiles in any horizontal plane, at 
a height y from the wall. To check the similarity of these profiles, 
u/U, was plotted against n, for different x stations. It was found 
that the velocity profiles in the transverse plane were all similar 
beyond a distance of about 6 to 8 times the nozzle height. Further, 
the distributions in different planes parallel to the wall behave 
alike and are similar. A typical non-dimensional distribution is 
given in Figure 4.5. It can be seen that the data can very well be 
described by Goertler's curve for a free circular jet (Abramovich, 1963). 
Figure 4.5 also shows that for planes close to the wall, the scatter 
is small and as we go higher, the scatter becomes appreciable. This 
again is attributed to the inaccuracies in measurement of low velo- 
cities. The data for all the shapes investigated behave very much 
alike, which is rather remarkable, and curves for the various cases 
studied are presented in Appendix C. The data are also compared with 
Tollmien's solution (Abramovich, 1963) for a free axially symmetric jet 
and an exponential curve. From Figure 4.5, one can see the Goertler's 
curve to be in closest agreement with the present experiments and 


therefore the same curve would be adopted in all future considerations. 


et | eae? 

ifs v0% bre Thame cag whase2 etnt <x soit 
sitions ait 19v0 ciotastettee zt avo We fs 

-ol sv muint xs sit -29269 ils nt tert nade at 37 pret as 2M 


-Fisw and mort \d zoni ar.0 sell somes eh + 


ie ; 
ts .snsia fsinostyon yas nt — 3 foof ay ont pdb heat 
.2otttovg szent to vtinelimta sit doors oT . Thaw Sls mort y idpted 5 
bruot 26w tl .enottese x Snovsttth w0oy .n Jenhages bajtolg 26w aie 
xeltmre [lo sY9w snsigq Ssetevensts ait nf esfftosq Wiroofev sds Srl 
Aedtwua .idpisd 6fsson sat e5mtt 8.03 3 Suods to sonstetb ¢ baoyed 
oyeod Tisw oft ot [affsxsq e5nslq tnsysttib ar anoriudiatetb ont 
2t mottudiatetb Isnotansmtb-non Tsotqys A aslimte o16 bas odtts 
ad [aw yrav aso sisb sat tent nos2sd fbo-ST .e.A ayer at nevip 
.(€d@f .fofvomeydA) Jo, xsiusvio seit & Tot svtv> 2'ysfivs0a yd badtypesb 
1sitso2 ot . [law ont oF s20fo eonsiq yot tent ewore oefe 2.8 ewer 
2tiT .sfdsioevggs 2amosed tetss2e ait .1Siptd op Sw es bane fieme et 


~ofav wol to tnsmeiwessm nf 2eafoswwooant ont of botudinsts at nfses 


\ 


floum, yisv evened botsptdeevnt 2oqsde ons ffs yot séeb oat. .2attto 
29260 evotisy ont yot eaviwo bas ,sldsdyemay verses 2f da ta atts 
ditw bsisqmoo o2ls ove etsb 9dT .9 xtbnoqqA nt bstnsesigq ie botbute 
tot aivismmye yfisixs so1t 6 sot (£00! .totvomsydA) nottyloe 2'ngiml [oT 
2'yot21809° sri: 952. ngs ano <2)P owner? mort | Jenaus Fettnonogxe ns bas 
bas ednant+agxs sne2e1g eda dttw ‘tngmeatps teazols at od of svi. 
-2notdevebtenos swaut Me at batqobs sd biuow sviuo singe ait aie 


74 


RUN No. 8 Xe x! 


i? CIRCULAR pNOZZLE Mies. coulg 753 
ara, © 12.0 0 38.7 
"y wo 4 173 a 440 
1.0 from ah v 22.6 v 49.4 
: ; * 28.0 a 
a oe * 60.0 
aa 


—— GOERTLER’S SOLN. 
a —-:— TOLLMIEN’S SOLN. 
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Drawing isovels in a future section would need the coordinates 
of the classical wall jet and the Goertler's curve, which are there- 
fore given in Table IV-2 for ready reference. 

The irregularities in the transverse velocity profiles, re- 
ported by Sforza and Herbst (1967) for slender orifices could not be 
detected for these tests: The reason could be that in the present 
tests the aspect ratios were not too large so that the velocity dip 
caused by the induced velocities of the surrounding vortex ring would 


be small and undetected by the present apparatus. 


4.5.3 Decay of Maximum Velocity 


Before attempting to find the maximum velocity at any section 
which incidentally is also the velocity scale, an important aspect - 
the virtual origin should be considered. The virtual origin is a 
ficticious point at which the velocity of flow approaches infinity in 
order to maintain a finite momentum flux from a point source. To 
establish this hypothetical point, Ta is plotted against x'/h (or 
x'/d) as shown in Figure 4.6, where x' is the distance measured from 
the nozzle section. For test runs in series I, measurements were 
carried to limited distances of x'/h < 64 and the results indicate that 
the virtual origin could lie anywhere between 2 to 8 times the nozzle 
height downstream of the nozzle. As the velocity profiles should in- 
deed approach similarity at very large distances from the nozzle, more 
extensive measurements were carried out in series II tests conducted 
inthe wide flume for distances of x°/h < 107... 1his series Of tests 


also elimated the side wall effects inherent in narrow flumes. The 
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TABLE IV-2 
COORDINATES OF THE SIMILARITY PROFILES 


Circular Free Jet 


Classical Wall Jet (Goertler) 


n (n) n f(n) 
O02 0.79 0 1.000 
0205 0.90 Cato 0.990 
0.08 ‘Onc }e) 0.20 0.964 
O10 0.98 0.30 0.930 
(Oa ie) HeOC 0.40 05.875 
0.20 0.995 0.50 0.820 
0.30 0.95 0.60 0.754 
0.40 0.890 0.70 0.692 
Oe50 Owes 0.80 0.626 
0.60 0.760 0.90 0.565 
0.70 0.690 1.00 0.500 
0.80 02625 T.0 0.450 
0.90 0.560 Tearg0) 0.400 
1.00 0.500 1.40 0.310 
Waite 0.450 1.60 0.240 
120 0.405 1.80 O2187 
Nasl@) 0.360 2.00 0.145 
1.40 0.310 2.20 021-55 
aR 8) 0.260 2.40 0.090 
1260 0.220 

P79 0.170 
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FIGURE 4.6 VIRTUAL ORIGIN (BASED ON VELOCITY SCALE) 
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results of the series II tests are also shown in Figure 4.6 and it 
can be seen that the virtual origin based on the far field data lies 
between 10 to 12 times the nozzle height downstream of the nozzle. 
Further, Figure 4.6 reveals that in the fully developed zone UpoX 
where x is the distance from the virtual origin. This relationship 
is in agreement with the one predicted from dimensional considera- 
tions and similarity analysis (Chapter II). 

Often investigators do not take into account the virtual 
origin, which might be justified to the extent that it avoids intro- 
duction of another unknown variable - the location of the virtual 
origin. However, it would not be right to compare the experimental 
results with any theoretical law, unless the distances are measured 
from the virtual origin. Sforza and Herbst (1967) presented their 
data making no provision whatsoever for virtual origin and found that 


Alle 


u__a(x') for all conventional orifices tested. The present results 


mo 

when plotted in a form similar to that of Sforza and Herbst, indicate 

that the exponent would lie between -1.14 and -1.20. The fact that 

the virtual origin is not located at the nozzle section is an in- 

herent property observed in most jets of finite size. As its location 

can vary over a wide range, the exponent of the velocity decay law too 
would vary accordingly. This seems to explain the reason why different 
exponents are reported for different investigations. The location of 

the virtual origin would depend on the nozzle contraction and to a greater 


degree on the level of turbulence in the approaching flow. Information 


on this aspect of the wall jet seems to be nonexistent, although, some 
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experimental results for plane free jets are available (Flora and 
Goldschmidt, 1969). 

From a practical point of view, it would be very useful to 
Show the variation of U0 with distance (x') measured from the nozzle 
section. In an attempt to establish a relationship which would in 
general be valid for any shape, Uno! Uo is plotted against x'/r', where 
r' is the hydraulic radius of the wall jet, defined as the ratio of 
the nozzle area to the perimeter of the fluid boundary at the nozzle 
(Subramanya, 1967). The results are plotted in Figure 4.7 and it can 
be seen that some noticeable shape effect still exists. An examination 
of Eq. 2.56 suggests that if the streamwise distances are non-dimension- 
alised by a length parameter Y A, where A is the area of cross-section 
of the nozzle, anunique relationship independent of the outlet shape 
should result. The experimental data are replotted in Figure 4.8 as 
Uno! Yo versus x'/Y A. and it can be seen that a single curve corre- 
lates the data extremely well. The curve is described by the follow- 


ing set of empirical relationships. 
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4.5.4 Length Scale b 
eer ees 


The variation of the length scale is shown in Figure 4.9. For 
both series of tests, bY grows linearily with the streamwise distance 
and is represented by a single relationship regardless of the nozzle 
Shape. The virtual origin for by would be an imaginary point where 
it tends to be zero, and in general it would be different from the 
velocity virtual origin discussed earlier. The results indicate that 
the virtual origin for the normal length scale is located behind the 
outlet at a distance of about 20 times the height of the nozzle. The 


data shown in Figure 4.9 can be represented by, 


b ; 
— = 0.90 + 0.045 “ (4.2) 


This rate of growth is comparable with that of Newman et al. (1971), 
found for a circular wall jet with air as the test media. From Eq. 4.2, 
the rate of growth normal to the wall is seen to be almost half that of 
a free jet and is because of the boundary restraint on the growth of 


large size eddies normal to the wall. 


4.5.5 Length Scale b, 


From the velocity profiles taken in any particular plane 
parallel to the bed, the length scale D3 was determined and its vari- 
ation is given in Figure 4.10. A total of four different planes were 
considered and Figure 4.10 indicates that b, grows linearly with x. 


This is particularly true for planes close to the wall. In addition, 
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the virtual origin for De is located in the downstream direction at a 
distance of about 6 times the nozzle width. The straight line vari- 


ation is represented by the equation, 


b ' 
Z 
= 0.20 = - 1.25 (4.3) 
where B = width of the nozzle. 


The above expression describes the experimental observations 
reasonably well for y/h = 0.5, whereas for higher values of y/h there 
appears to be some departure from this. Thus, b, depends on x as well 
as on y, although the y dependence is very weak. A high percentage 
of the total mass flow and momentum is contained within a small height 
from the bed, and therefore represents the zone of interest from a 
practical point of view. Hence,it is recommended that Eq. 4.3 be 
adopted to cover the whole range of y/h values. 

It is important to note that b. grows almost 4.5 times as 
fast as bY - Further illustration of this aspect is in Figure 4.11, 
which was obtained by injecting dye into the main stream. Newman et 
al. (1971) have also indicated this peculiarity for an air jet, the 
growth rates being comparable with the present values. For the series 
Il tests, the water is recirculated and eliminates any possible tempera- 
ture difference and consequent density differences. Therefore, the ob- 
served difference in growth rates of by and b, are attributed to the 
effect of the wall. The data from series I and II are in good agree- 


ment, giving an indirect indication that there was not any accountable 
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temperature difference for the first series of tests as well. 


4.5.6 Isovels 

In the earlier discussions, the form of the functional re- 
lations a and g, were found to agree closely with the classical wall 
jet and the Goertler's curve respectively. From an examination of 
Eq. 2.9, one can see that, u/U could be represented as a function of 
ny and ite and would be independent of x coordinate. From the co- 
ordinates of f, and gy given in Table IV-2, u/U 6 at a number of grid 
points were established. Isovels or lines of constant u/U ng were 
drawn by interpolation, and the resulting curves are given in Figure 
4.12. If a region near the wall is neglected, the isovels could be 
approximated by ellipses. This type of behaviour has been noticed 


earlier by Trentacoste and Sforza (1967) for three dimensional free 


jets. 


4.5.7 Centerline Boundary Shear 

The jet Reynolds number defined as Uy A/v is generally 
large for most cases encountered in practice. Hence, the right hand 
side of Eq. 2.63 can approximately be considered as a constant quantity. 
Combining Eq. 2.63 and Eq. 2.56, it can be shown that Col Ung /2 
is a constant for the fully developed flow region. The same conclusion 
can also be arrived at froma different type of analysis presented by 
Sforza and Herbst (1967). On the assumption that the bed shear dis- 


tribution in the transverse direction was similar, they derived an 


expression for the skin friction from the equations of motion as, 
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b q,-| 
se KK kx (4.4) 


where, K is an empirical constant, q, is the exponent in the X-Dy growth 
relationship. At large distances from the nozzle, qj approaches unity, 
and hence the skin friction coefficient in eq. 4.4 tends to attain a 
constant value. 

Figure 4.13 presents the experimental data for the various 
shapes studied. For x'/vY A greater than 50, for the range of Reynolds 


numbers investigated, the following is seen to be satisfied. 


T 
yr = 0.0065 (4.5) 
2 © Uno 
Hence,in this region 
a (4.6) 
om xe ; 


A more convenient type of plot is shown in Figure 4.14 in which 

tomlZ P te is plotted against x'// A , and a single curve is easily 
drawn through all data points. This plot is believed to be more use- 
ful as it avoids the virtual origin, and gives a clearer picture about 
the absolute values of shear stress close to the nozzle. In effect, 
the area close to the nozzle would need more elaborate consideration 


for protection against scour, since the intensity of shear stress drops 
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very quickly in a short distance. 


4.5.8 Off-Center Boundary Shear Distribution 


The distribution of the x component of boundary shear stress 
in the transverse direction is shown in Figure 4.3, and is for a circular 
nozzle 0.75" in diameter. The measurements were conducted during in- 
vestigation on swirling jets and detailed discussions are presented in 
Chapter VII. At this stage it would suffice to say that, the distri- 
bution is similar beyond a distance of 12 times the nozzle height 
and is very well described by Goertler's velocity distribution curve. 
The length scale for shear distribution b equals z at which the shear 


stress is Ton < and is found to vary linearly with x’. 


4.6 Summary 


From experimental observations on wall jets of various shapes, 
the velocity distribution normal to the wall and along the centerline 
was found to be similar and closely approximated by the classical wall 
jet curve. The velocity distribution in transverse directions was 
also found to be similar and agreed with Goertler's curve for a free 
circular jet. The maximum velocity at any section in the fully deve- 
loped zone was inversely proportional to the distance measured from 
the virtual origin. This is in accordance with the theoretical pre- 
dictions. However, to avoid the trouble of locating the virtual 
origin, an empirical curve was presented to predict uo This curve 
seems to be independent of the nozzle shape, if / A 1s selected as 


a characteristic parameter to non-dimensionalise the distance from 
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the nozzle. The length scales by and b, grow linearly, but at different 
rates. The virtual origin for the velocity scale and the two length 
scales were located at widely varying distances from the nozzle, which 
upsets any plan of comparison between theory and experiment. The 
distribution of boundary shear stress in the transverse direction was 
similar and followed the Goertler's velocity distribution curve. The 
centerline boundary shear coefficient defined by tom? © Ut was 

found to be a constant, and is around 0.0065. Convenient plots for 


estimating the centerline shear in terms of the nozzle efflux velocity 


Uy is presented, giving a clear picture about the shear stress decay. 
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CHAPTER V 
EFFECT OF SIDE WALLS 


5.1 Introduction 

The case of a three dimensional wall jet growing in an in- 
finitely wide channel is highly idealised and in most practical problems 
the channel is of finite width and the side walls would in some way or 
other affect the jet characteristics. This chapter presents an experi- 
mental study on the diffusion of a jet issuing from a circular outlet 
located symmetrically at the center of channels having finite widths. 
The bed and side walls were all smooth. 

Consider the case of a circular jet of uniform velocity Uy 
issuing tangentially on the bed of a channel of width W, with the tail- 
water depth sufficiently large so as to keep the outlet deeply sub- 
merged. On leaving the nozzle the jet expands, as in the case of a 
very wide channel, and depending on the width of channel the expanding 
jet would soon strike the side walls at some distance from the outlet. 
This would inhibit the transverse growth and consequently the decay of 
the maximum velocity. At the present time, there seems to be no in- 
formation available on this aspect. Abramovich (1963) mentions about 
the diffusion of a free jet in a finite space, which occurs in the 
case of an ejector. He finds that the process of equalization of 


flow parameters in a cylindrical mixing chamber occurs in such a manner 
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that the velocity field at each of its cross-sections appears as ie ane 
were the central part, bounded by the cylindrical walls of the chamber 
of the universal function, which expresses the dimensionless velocity 
field at the corresponding cross-section of a free jet. However, the 
flow field is significantly changed, as eddies form and result in 
backflow near the solid boundary. Apparently, it would be desirable 
to examine the present data on a wall jet, and see to what extent its 
behaviour would be comparable with the free jet. 

Some data on the diffusion of three dimensional jets in 
channels of finite width are available (Subramanya, 1967) for compari- 


son with the present results. 


5.2 Experiments 


A circular nozzle 0.375" in diameter was used in channels 
having different widths. The bed was a smooth galvenised iron sheet 
and the sides were fabricated from smooth painted plywood boards, the 
distance of which could be varied from 5.35 to 18.7 times the nozzle 
diameter. These channels of required width were stationed in flume 2 
and the bed elevation was adjusted to get a tangentially flowing jet. 
The velocity field was explored with the 0.05" diameter total head 
tube. The static pressures were measured by a Pitot-static tube lo- 
cated close to the total head probe. No pressure traverse was made 
normal to the wall and it was assumed that the pressure distribution 
is hydrostatic. The same total head tube was used to measure the 


boundary shear. 
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5.3 Results 
5.3.1 Decay of Maximum Velocity 

The maximum velocity at any section occurred on the center- 
line and its variation is shown in Figure 5.1] as WAN Fa Vse Xu) deo 
is seen that Vauee varies linearly with the distance, for x'/d less 
than about 35. Moreover, the efflux velocity being the same for all these 
runs, from Figure 5.1, it can be said that the velocity decay law is 
the same in this zone irrespective of the channel width. It is to be 
noted that this linear law near the outlet is different from the one 
derived for x'/d greater than 35, in case of a very wide channel. For 
x'/d greater than 35, the channel width appears as a governing para- 
meter and the velocity decay follows different laws. 

For better comparison, the data are plotted with Uno! Uo 
versus x'// A in Figure 5.2a. For values of x'/y A less than about 
45 (which also corresponds to x'/d = 35) the relationship is inde- 
pendent of the channel dimension. It is also interesting to note that 
the data of Subramanya (1967) on rectangular three dimensional wal] 
jets of different aspect ratios (in the range of 0.94 to 34/2 issuing 
into channels of finite but different widths, follow closely the mean 
curve of the present investigation (Figure 5.2b). His observations 
were limited to x'/v A < 27, on the basis of which he concluded that 
the relative width of the channel does not appreciably affect the 
velocity decay. This is in accordance with the present observations 
on a circular jet. However, at large distances from the nozzle the 


entire flow would be in the forward direction and tend to behave as 
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an open channel. This would mean that the maximum velocity at any 
section would tend to attain a constant value. This trend can be seen 
in Figure 5.2a, and is particularly clear for the narrowest channel 
studied. The higher values of Uno! Yo is also evident from the fact 
that the supply discharge remains same for all tests and therefore 


narrower channels will have a higher forward flow velocity. 


5.3.2 Length Scale b 
ein ETE) 


The length scale found from the velocity profiles taken along 
the centerline is given in Figure 5.3. Here again, the variation is 
linear for x'/d less than about 35 and is practically independent of 
the channel width. The specified limit for x'/d = 35 increases with 
an increase in the width, beyond which the growth rate is very fast 
and the change is abrupt. In the first instance, this was thought 
to be caused by the adverse pressure gradients. The measured water 
surface profiles are given in Figure 5.4, and as such the effect of 
the adverse gradient should not be felt for x'/d less than 60 to 70. 


No other suitable explanation could be offered for this peculiarity. 


5.3.3 Velocity Profile Normal to Bed 


The profiles taken along the centerline in a normal direction 
were checked for similarity by plotting the conventional quantities, 
Un! Uno versus y/by» in Figure 5.5. The scatter in the data was large 
and the velocity distributions in the narrow channels were not similar. 
As the channel gets wider, the distribution tends to follow the classical 


wall jet profile. This is seen from the tests on the channel having 
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FIGURE 5.4 WATER SURFACE PROFILE 
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W/d = 18.7, which is the maximum width in this series. It is hard to 
demarcate precisely the limiting value of W/d beyond which the velocity 
profiles could be considered as similar. However, for the tests con- 
ducted in flume 1 while investigating other aspects of the jet, the 
velocity profiles were found to be similar, the corresponding W/d value 
being around 33. Therefore, as a rough guide a limiting value of 30 is 
Suggested, beyond which the results presented for a jet in a wide 


channel (Chapter IV) could be applied successfully. 


5.3.4 Transverse Distribution 

Some velocity measurements were taken in the x-z plane lying 
close to the bed. At a short distance of 5.3d, the velocity profiles 
look similar to the one found for a wide channel, but soon the entire 
width is occupied by the forward flow and the half-width could not be 
estimated, as a result of near uniform distribution of velocity. 
Measurements at a distance of about 35d, for the 2" wide channel for 
example, show a transverse distribution that is uniform, and the same 
characteristic is observed for any section beyond this limit. In general 


all the widths tested behave in a similar fashion. 


Saoeomepedrollcar 

The bed shear stress along the center of the channel was 
measured by the Preston tube technique. From velocity measurements 
taken close to the wall, and the shear estimated by a Preston tube, a 
check was made on the validity of the wall-law. The logarithmic law 


was found to be reasonably satisfied for all the tests and therefore 
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the shear measurement technique is satisfactory for narrow channels as 
well. 

The variation of the local skin friction coefficient is shown 
in Figure 5.6. The results are compared with the mean curve for three 
dimensional jets on smooth beds in wide channels. The agreement looks 
fair for x'/vV A less than seventy five. Beyond this distance, Cr 
drops from the constant value of 0.0065 and it is suspected that the 
adverse pressure gradient acting in this zone causes this. 

Same experimental data are presented by defining a skin 
friction coefficient in terms of the nozzle efflux velocity, and this 
alternate plot is shown in Figure 5.7. For x'/v A greater than about 
40, the coefficient is much higher for the channels having finite 
width. This would mean that, under identical efflux conditions of the 
nozzle, the centerline bed shear is higher when discharging into finite 
width channels. As the channel gets narrower, the shear intensity 
increases. As an example, for the channel having W/d = 5.35, the 
intensity of bed shear at a reference section of x'/¥ A = 100 is 
around four times larger than the corresponding value for the wide 


channel case. 


5.4 Summary 


In a certain initial reach, the velocity decay and bY growth 
relations do not get affected by the channel dimensions. For narrow 
channels the maximum velocity along the centerline tends to attain a 


constant value, higher than what would be if it were to diffuse in an 
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FIGURE 5.6 SKIN FRICTION COEFFICIENT (ALONG THE CENT ER -PLANE) 
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infinitely wide channel. The bed shear for identical outlet conditions 
would be higher in the far field, for channels having finite width. 

As the channel gets narrower, the flow field is further complicated 

by the action of pressure gradients, and similarity of velocity 
profiles is not satisfied. A limit for W/d is roughly set at 30, 


beyond which most results of the wide channel investigation seem to 


apply. 
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CHAPTER VI 
EFFECT OF BOUNDARY ROUGHNESS 
ON CIRCULAR WALL JETS 


6.1 General 

Studies on diffusion of a circular jet impinging tangentially 
on a rough boundary, under zero pressure gradient, are presented in this 
section. The ambient fluid was stagnant and the outlet was deeply sub- 
merged so that the ambient could be treated as approximating a semi- 
infinite medium. It has been established in Chapter IV that on a 
smooth boundary all shapes of nozzle behave alike; moreover, the 
circular shape being more common in hydraulic engineering, investi- 
gations were carried out for this shape only. 

The other investigations on wall jets growing on rough 
boundaries appear to be those of Lam Lau (1963) and Rajaratnam (1965). 


These treated only the two dimensional case. 


6.2 Existing Work on Plane Jets 


Rajaratnam conducted experiments on plane turbulent wall 
jets under zero pressure gradient on rough boundaries and covered a 
range of K/h varying from 0.00455 to 0.126, h being the initial 
depth of wall jet. It was observed that in the fully developed zone, 
the velocity profile should be considered in two parts, viz, the 


inner or the boundary layer and the outer or free-mixing region. Al- 
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though a single length scale does not exist, each of the parts exhibit 
similarity if the appropriate length scales are considered individually. 
The variation of the two length scales and the common velocity scale 

Uno Were presented. A power law type of equation described the velocity 
distribution in the boundary layer, with its exponent varying from 1/2 
to 1/8. The decay of the maximum velocity was reported to be affected 
by the boundary roughness. Studies on momentum fall showed that for the 


K/h = 0.126 case, the flow could loose more than 50% of its efflux 


momentum by the time it travels a distance of 60h. 


6.3 Dimensional Considerations 

For flow over rough boundaries the loss in forward momentum 
is high and it is not possible anymore to assume that the efflux mo- 
mentum flux is preserved. This would mean that the efflux momentum 
flux can no longer be considered as a single variable, but the efflux 
velocity Uy and the area of the nozzle would appear as separate vari- 
ables, affecting the flow field. Hence, for turbulent flow on a rough 


boundary, the maximum velocity at any section may be expressed as 


Vier a F(Up »A.p XK.) (6.1) 
u K 

or Toe (one (6.2) 
0 YA VA 


For flow on a smooth boundary it has been established (Chapter II) 


that, 


he a 
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mo _ G 
li 3 Se (2.56) 


0 (x/v A) 
An examination of Eqs. 6.2 and 2.56 (together with the experimental 
results presented later) suggests that Eq. 6.2 could possibly be re- 


written as 


U C(K./v A) 
_=— (6.3) 
0 xX/V A 
and for a circular outlet 
u C(K_/d) 
Ae aes (6.4) 
U5 (x/d) 


As data on plane wall jets on rough boundaries are extensive, it would 
be desirable to write the corresponding expression for the plane jet 


as 


uU C(K./h) 
mo = —2__ (6.5) 


U5 (xh)o-? 

The data of Rajaratnam were analysed by plotting Uno! Yo versus x/h on 
a double-log plot. By drawing lines at a slope of 0.5, which would 
give best fit for the test data for any particular run, C was deter- 
mined for various roughness heights. The variation of C for a wide 


range of relative roughness is shown in Figure 6.1. There is no systematic 


variation although there seems to be lot of scatter and C remains un- 
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hy 
affected by the roughness. This is contrary to the belief that C should 
decrease with an increase in bed roughness. Observations on the plane 
jet indicate that an increase in the roughness displaces farther away 
from the boundary the point of maximum velocity and the distribution 
in the boundary layer is less steep. This distribution of velocity, 
added with the thicker boundary layer, account for the high rate of 
momentum loss which is almost entirely confined to the wall region. 
The wall does not affect the free-mixing zone and the maximum velocity 


occurring at any section to any significant extent. 


6.4 Present Investigation 


By varying the boundary roughness systematically, design 
curves for predicting the gross behaviour could be prepared. However, 
the discussions on the plane jet (Section 6.3) dictate that the present 
investigation should aim at finding out whether the results for the 
smooth boundary are applicable as an approximation to the rough bo- 
undaries. If not,the existing order of differences should be found. 
Other details about dividing the flow field into the inner and outer 
region, the velocity distribution, boundary layer growth along the 
centerline and the decay of maximum shear were also studied for a 


circular wall jet. 


figt Experiments 
Experiments were carried out in flume 2. The rough boundary 


was formed by gluing a continuous strip of #36 wet-or-dry aluminium 


oxide cloth (described in Section 3.4) on the smooth aluminium plate. 
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The nozzle was 0.375" in diameter, and the jet location was adjusted 
till visually it could be judged to be issuing tangential to the rough 
boundary. For this test run, the relative roughness K./d = 0.195. In 

a major portion of the flow, the velocity field was explored by a total 
head tube 0.085" inside diameter and 0.11" outside diameter. However, 
detailed mean velocity measurements in the boundary layer along the 
centerline and the velocity measurements close to the nozzle were taken 
with a 0.05" total head tube. The static pressure was measured by means 


of a side-wall tapping. 


6.6 Results 
6.6.1 Similarity of Velocity Profile 

Non-dimensional plots Ue! Uno VS. y/by are shown in Figure 
6.2 and indicate similarity of the profiles, for x'/d in the range of 
15 to 50. Beyond the upper limit, the scatter increases and seems to 
be the result of inaccurate mean velocity recorded in this zone of 
high turbulence intensities. Comparison of the experimental similarity 
curve with the classical wall jet shows that the rough wall shifts 
outward the point of maximum velocity. Also, use of a single length 
scale by for the entire thickness gives wide scatter for the points in 
the boundary layer, and it should be checked as to whether division into 
two regions having different length scales would improve the velocity 


distribution for a three dimensional jet. 
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6.6.2 Free Mixing Region 

The shear layer lying above the boundary layer of thickness 
6 is referred to as the free mixing region. The distances and the 
velocity scale for this zone were measured from the outer edge of the 
boundary layer. Hence, the non-dimensional velocity distribution is 
given by f(n') (= Un! Umo) VSenm (= y-6/by~6) in Figure 6.3. In the 
Same figure, the data on the smooth boundary is plotted for compari- 
son. In both cases, the distribution agrees closely with the classical 
wall jet curve. Thus the free-mixing region is unaffected by the 


particulars of the boundary roughness. This is identical with the 


results of Rajaratnam (1965) for plane wall jets. 


6.6.3 Inner Region 

From a knowledge on the velocity distribution in the familiar 
boundary layer and pipe flows, the present data are plotted as U Ung 
versus y/6 in Figure 6.4. For small distances from the nozzle, K./6 
would be large compared to the values for large distances and therefore 
a Single curve is not to be expected. Data indicate a two fold vari- 
ation in the value of K./6. However, Figure 6.4 shows that excluding 


the first section, the rest of the data can be described by a simple 


power law distribution, 


A value for n = 8 in the above expression gives a reasonable fit to 


the data. Rajaratnam found that for a plane jet, the value of n can 
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FIGURE 6.4 VELOCITY DISTRIBUTION IN THE BOUNDARY LAYER 
(ROUGH WALL) 
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be between 2 and 8 and the lower value occurs for the higher roughness. 
For a three dimensional jet, the growth of the boundary layer normal 
to the bed is limited by the flow spreading in a transverse direction, 
and therefore results in a high value for n even for case of a very 


rough wall. 


6.6.4 Velocity Distribution in x-z Plane 

Velocity distribution in the transverse plane is shown in 
Figure 6.5 and represents the conventional quantities u/u., VS. z/b_. 
The distribution is found to be similar and for planes lying close to 
the boundary, the agreement is good. Goertler's free circular jet 
curve is compared with the data and can be seen to agree closely 
with the experimental results. Thus the boundary roughness does not 


affect the nature of this distribution to any noticeable extent. In 


planes lying at higher elevation from the bed, the data scatter more. 


6.6.5 Maximum Velocity Decay 


The maximum velocity Ur along the centerline is a scale 
common to both the inner and the free mixing regions. The observations 
on the maximum velocity are given in Figures6.6. Plot of a. versus 
x'/d in Figure 6.6a for the rough and the smooth boundary, gives linear 
relationships. The two tests were conducted for identical efflux 
velocities and therefore equal slope of the two linear variations 
would mean that the rate of decay is the same, regardless of the 
boundary roughness. On a rough boundary the process of mixing is en- 


hanced and results in a shift in the virtual origin in the up-stream 
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‘direction. Thus,at a section equidistant from two identical outlets, 
the velocity for the rough boundary case would be somewhat less than 
the smooth wall case. However, from the following discussion, it 
seems the differences in the two velocities are small. 

As an alternative, the same data are shown in Figure 6.6b, 
as Uno! Yo vs. x'// A. The mean curve derived from tests on all nozzle 
shapes is seen to describe the rough wall data extremely well. This 
observation, combined with the analysis of data (Section 6.3) for 
plane jets on rough walls, suggests that the maximum velocity decay 
law is the same for all kinds of boundaries. Hence, the relationships 
derived from extensive investigation on three dimensional wall jets 
on smooth walls can satisfactorily be used in practice, unless the 


relative roughness is excessively large. 


6.6.6 Length Scales 


As the centerline velocity distribution normal to the wall 
is considered in two parts, two length scales need to be specified. 
In addition, there would be a third length scale for the transverse 
distribution. 

The length scale for the free mixing region is (b-6), and 
its variation with distance x' is given in Figure 6.7a. The variation 
can be found to be linear and follows closely the data on smooth boundaries. 
Beyond a distance of x'/d > 50, the length scale for a rough wall grows 
fast, and no explanation could be found. A general relationship which 
would be valid for all types of boundary roughness is given by 


(b,-6) 


—Y,— = 0.042 (4) + 0.5 (6.7) 
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Variation of the boundary layer thickness 6, which is the 
length scale for the inner region, is given in Figure 6.7b. For the 
smooth as well as the rough boundary the variation of 6 with x' is 
linear, the values being higher for the rough boundary. It is to be 
pointed out that in all tests it was impossible to locate the outlet 
exactly tangent to the boundary and a small step is provided at the 
outlet section. As this step size is arbitrary, and is likely to be 
different in the two sets of experiments, no significance is given to 
the boundary layer thickness measured close to the outlet section. 


The growth of the boundary layer for a circular outlet is given as 


$= 0.01 (4) +k (6.8) 


where K is a constant and depends on the boundary roughness and the 
nozzle contraction details. 

Growth of b. in three different planes was considered and 
its variation is given in Figure 6.8. Up to x'/d = 30, the growth 
rate is very much like the smooth boundary case, but beyond this 
distance the growth rate of b, is faster. Even the data on a plane 
close to the rough boundary does not exhibit linear variation. There- 


fore, no simple expression is suggested and b. should be read from the 


plot. 
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FIGURE 6.7 GROWTH OF LENGTH SCALES FOR FREE MIXING AND INNER REGION 
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G..0°. 7 Boundary Shear 


Beyond an initial reach of 2.5" from the outlet, the center- 
line bea shear was computed from Preston tube readings. A larger tube 
of 0.00917' outer diameter was used to minimize the error in location 
of the zero datum. From Eq. 3.11, the centerline bed shear was found, 
ana the values seemed very high and unrealistic. In the range of 
distances considered, the tube occupies 30 to 75% of the boundary layer 
thickness, so that, the displacement of the effective center would be 
large compared to the boundary layer thickness. This would mean that 
in all cases, tne total head tube records a much higher pressure read- 
ing than actual. 

The velocity profile in the boundary layer was then measured 
with a smaller total head tube of 0.0041' outside diameter. The data 
plotted as straight lines on a semi-log paper and the bed shear was 
computed from Eq. 3.13. When the tube is resting on the boundary, the 
displacement effects would be a maximum, and therefore, the first point 
should be avoided as far as possible in fitting a straight line. The 
boundary shear found by this method was found to be low, and in many 
regions the local skin friction coefficient was found to be somewhat 
smaller than the one for a smooth boundary, which is rather baffling. 

Until more reliable results are obtained for boundary layers 
of larger thicknesses (using bigger nozzles), no comment could be 


made on the bed shear distribution on a rough boundary. 
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6.7 Summary 

The centerline velocity profile normal to the boundary can 
be split into two regions - the free mixing zone, and the inner zone, 
which exhibit similarity. The distribution in the free mixing region 
closely follows the classical wall jet profile, regardless of the 
boundary roughness. A power law type of distribution is found to be 
satisfactory for the inner region. The length scale for the free mixing 
region and the boundary layer thickness, vary linearly with distance x. 


Pray) and the boundary roughness does not 


The maximum velocity Uno & * 
affect the x'/vV A vs. Uno! Yo curve to any significant extent, so that, 
for preliminary estimates, the curve for the three dimensional wal] 
jet on a smooth boundary would be used. 

The velocity distribution in the transverse direction can be 
approximated by Goertler's curve, noting that the growth of b, can no 
longer be represented by a linear relationship. Both the Preston tube 


technique and velocity profile method, proved unsuccessful in giving a 


true picture of the centerline bed shear distribution. 
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CHAPTER VII 
WALL JETS WITH SWIRL 


7.1 Introduction 

Giving a certain amount of swirl or rotation, before a jet 
comes out of a circular nozzle, helps to spread the flow more rapidly 
and therefore the maximum velocity decays faster. Considerable amount 
of work has been done on the effect of swirl on free circular jets and 
these swirling jets find application in combustion chambers and 
sprayers. Experimental observations on free swirling jets have been 
made by Rose (1962), Gore and Ranz (1964), Kerr and Fraser (1965), 
Chigier and Chervinsky (1967), Pratte and Keffer (1969) and others. 
In the fully developed flow region the distribution of axial and 
tangential velocity components were found to be similar, with the 
length scales varying linearly with distance x. The pressure plus 
momentum (P+M) and the angular momentum (T) were found to be preserved 


and were combined to form a dimensionless parameter, 


atl. 
where S = swirl number, and 
Na) = radius of the nozzle 


From a similarity analysis it is possible to show that for a free 
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swirling jet, when the maximum axial velocity Hi is far larger than 
the maximum tangential velocity cies at the nozzle, the following re- 


lations are satisfied. 


bax 


(de2) 


where b is a length scale. 


From the experiments of Chigier and Chervinsky, even for 
the hignest value of S(= 0.64) the above expressions were found to be 
valid. 

There appears that no information is available on the effect of 
swirl on a circular wall jet. In an attempt to find how effectively 
one can spread the flow from a circular outlet resting on a smooth 
boundary, some preliminary experiments were conducted. At this stage 
any reasonable theoretical approach is not possible because of lack 
of symmetry in the flow field. In this chapter the forward velocity 
and the bed shear distribution for S = 0, 0.141, and 0.265 were measured 


in detail. 
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7.2 Experiments 


Experiments were conducted in flume 2 for a circular nozzle 
0.75" in diameter. Selecting proper axial and tangential discharges 
through the swirl generator, the swirl number could be suitably adjusted 
in the range of 0 to 0.265. The velocity field was explored by a 
calibrated five hole probe 0.17" thick, and the bed shear distribution 
was explored by a yaw probe. The tail water depth was of the order of 


28 times the nozzle diameter and resulted in a level surface. 


7.3 Results and Analysis 
7.3.1 Velocity Distribution Normal to the Wall 


Distribution of the axial velocity component taken normal to 
the wall are given in the non-dimensional form ny vs Un! Umno (Figure 
7.1). Observations were limited to distances of 40 and 26 times the nozzle 
diameter, for S = 0.141 and 0.265 respectively, since beyond this 
distance the velocity components were smal] and would have introduced 
large errors in measurement. For both the experiments, the distribu- 
tions were found to be similar and as in the case of non-swirling jets, 
the distribution is closely represented by the classical wall jet curve. 
For swirling jets, similarity was found to be satisfied beyond a short 
distance from the outlet. For example, for S = 0.265, velocity pro- 
files were similar beyond a distance of 6d from the nozzle. These 


observations were in general similar with those on free swirling jets. 


7.3.2 Velocity Distribution in the Transverse Direction 
Distributions of axial velocity components are indicated in 
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Figure 7.2. For y/d = 0.139 and 0.512 the profiles look similar and 
are very well described by Goertler's curve for free circular jets. 
Measurements in higher planes showed a shift in the location of the 
maximum velocity from the centerline. This shift was of varying degrees 
depending on the distance of the section from the nozzle. Also, many 
of the profiles looked unsymmetrical and irregular. Measurements at 
large distances from the nozzle proved unsuccessful because the velocity 
components were small in the higher planes. As a result, the investi- 
gation had to be limited to locations close to the bed, and the dis- 
tributions were found to be similar beyond a short distance of 4 to 6d 


from the nozzle. 


j-o2d sDISTYIDUTION Of Voandiw 

As a typical example, the distribution of v and w components 
on the centerline normal to the wall are given in Figure 7.3 for 
S = 0.141. The w component decays very fast and the distribution 
does not look similar. The v component was also unsymmetric and the 
present data prevent us from drawing any definite conclusion. The 
same is true for the static pressure distribution (Figure 7.4). Be- 
yond a short initial reach, all the above mentioned quantities are of 
very small order and their measurement with the present technique is 


likely to be not very reliable. 


7.3.4 Decay of Maximum Velocity 


Maximuin velocities measured along the centerline of the jet 
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FIGURE 7.2 NON-DIMENSIONAL VELOCITY DISTRIBUTION IN TRANSVERSE PLANES 
OF SWIRLING JETS 
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FIGURE 7.3 v AND w VELOCITY COMPONENTS - A TYPICAL PLOT 
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were non-dimensionalised with respect to the average efflux velocity 
at the outlet, and the variation is given in Figure 7.5. The data for 
zero swirl on the 0.375" and 0.75" nozzle were in close agreement. 
Velocity measurements at the 0.75" nozzle efflux section showed that 
the distribution was not uniform, but resembles the turbulent pipe flow 
profile. This is perhaps because of the straight length of 1" dia- 
meter pipe provided upstream of the test nozzle and the insufficient 
contraction ratio at the nozzle entry. First, the maximum velocity 
at the efflux section was selected as the reference velocity Uy and 
the data showed no agreement with the results from the 0.375" dia- 
meter nozzle which had an uniform velocity distribution at the out- 
let. Hence, from the measured discharge rate the average velocity 
was computed for the 0.75" nozzle and was taken to be the representative 
velocity Uy: This seems justified in the light of agreement between 
the results for different size nozzles. 

Even for jets with swirl, Up was the average forward velocity 


and the decay of the maximum velocity Un is shown in Figure 7.5. The 


) 
data indicate u. a 1 for this range of swir] numbers and that the 


mo 
decay is faster for the swirling jets. One can clearly see that a 
small swirl is very effective and for higher swirls the effects do not 
increase proportionately. The difference in the decay rates for S = 0 
and 0.141 is appreciable, whereas, the difference for S = 0.141 and 0.265 
is very small. Another important feature of swirling jet is that the 


flow becomes similar at relatively short distances from the nozzle, 


and the virtual origin shifts close to the outlet section. As a result, 
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FIGURE 7.5 MAXIMUM VELOCITY DECAY FOR SWIRLING JETS 
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the overall effect of swirl is to have at identical sections velocities 


far less than those for non-swirling jets. 


7.3.5 Growth of Length Scales b. and b, 
hee ee ee oe ee) AS ee 


For a short initial length of 3d, by does not grow at all and 
then falls for a swirling jet. From Figure 7.6, this dip can be seen 
to increase with higher swirls. At around 8d, the jet starts to grow 
and beyond about 20d the growth becomes linear and follows the data 
for non-swirling jets. Hence, in the far field (x'/d > 30) the re- 
lationship developed for three dimensional wall jets (in Chapter IV) 
may be used as a approximation for jets with moderate swirls. 

For comparison, growth of b, close to the wall is shown in 
Figure 7.7, for different swirl numbers. Here again, the expansion 
rate is more for swirling jets which increases with swirl number. 
Also, jets with swirl tend to expand right from the outlet section 
causing the virtual origin to lie very close to the nozzle (Figure 
7.7). For the range of swirl numbers investigated, the growth of b, 
close to the wall is linear, which is typical of the behaviour of b, 
in higher planes as well. In the above discussion no attempt was made 
to give empirical expressions for by and b> since we have results 


for only two non-zero values of S. 


7.3.6 Bed Shear Distribution 
Transverse distribution of the x component of the bed shear 
was measured at a number of x stations and it was observed that the 


maximum value occurs almost on the centerline and the distribution was 
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more or less symmetrical. The maximum shear stress eae UE the center- 
line was taken as the scale for shear stress and b. was taken as a 
length scale, being equal to z at which Ce Tom/*: The non-dimensional 
Shear distribution is given by TaslTeea vs. ele z/b_) as indicated in 
Figure 7.8. For all cases, the shear distribution is similar in the 
fully developed zone and is very well described by Goertler's velocity 
distribution curve. For higher swirl numbers the distribution be- 


comes similar very quickly at a distance of 3d from the nozzle. 


A gaussian profile given by 


16 
-2_ = exp -{0.693 n,* } (7.3) 
om z 
is also indicated in Figure 7.8 and this empirical relationship does 


not give any improved fit to the data. 


7.3.7 Length Scale b 


From the present investigations, the shear length scale is 
found to vary linearly with distance x. Also the virtual origin is 
located close to the nozzle and the growth rate increases with the 


degree of swirl (Figure 7.9). 


7.3.8 Alternative form of Bed Shear Distribution 
Since the virtual origin for b. lies close to the nozzle, 
jt was thought that use of the shear length scale could be avoided 


if z is non-dimensionalised by the distance x'. The alternative way 
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of representing the bed shear distribution is given in Figure 7.10. 

In case of S = 0 and 0.141, the distribution for x'/d < 10 follow 
different curves and beyond this limit, all the test points tend to 

fall into a single curve. For the highest swirl (S = 0.265), the 
distribution becomes similar even at x'/d = 3.0. Experimental curves 
are drawn through the points lying on profiles that exhibit similarity 
and are then assembled together in Figure 7.11. This indicates that 

the shape of the similarity distribution changes with S and the lateral 
extent in which the bed shear acts increases greatly with an increase in 


Swirl. 


Jvaeo) Decay or Center] ines bed Shear Stress 


The centerline bed shear stress represents the maximum at 
any section. This is non-dimensionalised by the mean efflux velocity 
and the streamwise variation is given in Figure 7.12. For the case of 
zero swirl with the 0.75" diameter nozzle the centerline bed shear 
was determined by a total head tube as well as with a yaw probe, and 
the results are shown in Figure 7.12. The agreement of the data 
from the two methods was close and both follow the mean curve derived 
earlier for three dimensional wall jets of various shape (Chapter IV). 
This gives confidence in the yaw probe technique for shear measure- 
ments. For jets with swirl con falls very rapidly right from the out- 
let section. The rapid fall in the intensity is expected as the la- 
teral extent over which the bed shear acts increases very much with 


swirl and would therefore give a far smaller peak value on the center- 
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FIGURE 7.10 NON-DIMENSIONAL BED SHEAR DISTRIBUTION 
(WITH x' AS THE LENGTH SCALE) 
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FIGURE 7.12 DECAY OF CENTERLINE SHEAR STRESS FOR SWIRLING JETS 
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line. Here again, a weak swirl is found to be most effective and the 
effects are not enhanced very much for stronger swirls. The maximum 
shear stress at corresponding sections of swirling jets is 2 to 2.5 


times less than that for a non-swirling jet. 


7.3.10 Transverse Bed Shear Component 

In the discussions so far, only the x component of bed shear 
vector has been considered. However, some knowledge is needed about 
the z component if one is to have a complete picture of the shear 
distribution below outlets. Of all the alternative ways of repre- 
sentation, it was felt that an indication as to how the yaw angle would 
vary with z would be the simplest to get a rough estimate of the trans- 
verse shear component, once the x component u is known. Hence, for 
the various swirls investigated, the data for the yaw angle 6 vs. 2/x' 
were plotted for a number of x stations (Figure 7.13). As the swirl 
component at sections close to the nozzle were strong but unsymmetric, 
because of the wall, 9 does not behave systematically. For distances 
beyond 8 to 10d, the flow looks symmetric about the x-axis and 6 


increases linearly with z/x'. Further, the data from various sections 


collapse into a single curve. 


7.4 Summary 
The velocity distribution along the centerline of the jet 


normal to the wall was found to be similar. The distribution of axial 
components of velocity and bed shear in the transverse direction were 


also found to be similar and closely followed Goertler's curve. Even 
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for the highest value of S = 0.265, all the length scales and the maxi- 
mum velocity decay followed a linear variation with x. The growth rate 

of the jet and the velocity decay were faster for swirling jets. Moderate 
swirls were found to be very effective in reducing the centerline bed 


shear stress. 
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CHAPTER VIII 
CONCLUSIONS AND RECOMMENDATIONS 


8.1 Summary and Conclusions 


In hydraulic engineering, flows from deeply submerged circular 
and rectangular outlets are often encountered. A turbulent wall jet 
issuing from a circular, elliptic, triangular or rectangular outlet of 
finite aspect ratio is three dimensional in nature. As an initial 
step towards understanding the mechanics of flow from submerged out- 
lets, a simple free jet was considered for analysis. Based on the 
observation that the forward-momentum flux of the wall jet was not 
reduced to any appreciable degree, the free jet analysis was extended 
to three-dimensional wall jets. Extensive experimental investigations 
were then conducted on wall jets issuing into a wide channel from 
nozzles of various shapes. Aspects of practical interest, such as the 
effect of the channel side walls, the boundary roughness, and swirl 
were experimentally investigated on circular wall jets. A summary of 
the investigation on each aspect is included in the respective chapters. 
The significant conclusions drawn from the present investigation are 
collected and are presented below. 

From a similarity analysis (Chapter II) on the boundary- 
layer type of equations, the maximum velocity decay and the growth 
of the length scales for a three dimensional wall jet were given by 


the following expressions. 
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Dimensional considerations also indicated the same type of variations. 
Moreover, dimensional considerations indicated that Y A , the square 
root of the cross-sectional area of the outlet (Eq. 2.63), should be 
taken as a characteristic length if the shape effects of the nozzle 
are to be eliminated. An expression for the centerline bed shear stress 
variation was derived in terms of the characteristic parameter /Y A and 
the Reynolds number of the jet at the efflux section. 

Experimental observations (Chapter IV) on wall jets of 
various shapes showed that in the fully developed flow (x'/h > 15), the 
velocity profiles satisfied similarity and the classical wall jet curve 
served as a good approximation for describing the data in the center- 
plane. The distribution in the transverse direction of the fully- 
developed flow were also found to be similar, following closely Goertler's 
free circular jet curve. The theoretical relationships given in Eq. 
8.1 were in general found to be satisfied. The length scale b, grows 
4 to 5 times faster than by: The virtual origins for the velocity 
scale and the two length scales were located at widely varying dis- 
tances from the nozzle, suggesting that the flow was not truly self- 


preserving. Use of Y A to nondimensionalise x', as suggested by dimen- 


(1.8) 
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sional considerations, describes the maximum velocity decay relation- 
ship by a unified curve regardless of the nozzle shape. The boundary 
shear stress distributions in the transverse directions were similar 
and followed Goertler's velocity distribution curve. The centerline 
boundary shear stress coefficient cel = wad F ung”) equalled 0.0065 
in the far field and was found to be a constant. 

It was observed in Chapter V that when a circular wall jet 
diffuses in channels of finite width, in a certain initial region 
(x'/d < 35, for the range of widths considered in the present investi- 
gation) the velocity decay and growth of by were not affected by the 
side walls of the channel. However, the similarity of the centerline 
velocity profiles was not satisfied. The maximum forward velocity and 
the centerline bed shear stress in narrow channels were higher compared 
to the values for a infinitely wide case. A lower limit for W/d of 30 
is tentatively suggested, beyond which most results derived from investi- 
gation on a three dimensional wall jet discharging into a wide channel 
seem to apply. 

For flow over a rough boundary (Chapter VI), the center- 
line velocity profile normal to the wall was considered in two parts - 
the free-mixing zone and the inner zone. The velocity distribution 
in the free-mixing region was similar and was very well described by 
the classical wall jet profile regardless of the boundary roughness. 
The distribution in the inner zone was also similar and a power-law 
type of representation was satisfactory for representing the data. 


The growth of length scales for the free-mixing region was unaffected 
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by the boundary roughness and the streamwise variation was linear. 


4 and the boundary roughness did 


The velocity scale Ue varied as x 
not seem to affect the relationship x'/V A vs. Uno! Yo to any significant 
extent. This finding was also supported by the analysis of the avail- 
able data for plane wall jets on rough boundaries (Section 6.3). Hence, 
for preliminary estimates, the curve for the three-dimensional wal] 
jet on a smooth boundary is recommended for use in engineering practice 
so long as the roughness of the bed is not excessive. The growth of 
the length scale b, on a rough boundary was not linear, although the 
distribution in the transverse direction was similar and was given by 
Goertler's curve. 

Experiments on circular swirling wall jets (Chapter VII) 
with moderate swirl showed the centerline velocity distribution normal 
to the wall to be similar and could be described by the classical wall 
jet curve. For moderate swirls, all the length scales, including the 
one for bed shear distribution, followed a linear variation with x. 


steel The transverse growth rate 


Also, the maximum velocity Unio % X 
of the jet and the decay of the maximum velocity were faster for swirling 
jets, and in particular moderate swirls proved to be very effective. 
Swirls were also found to be very effective in reducing the centerline 


bed shear stress. 


8.2 Recommendations for Further Study 


For better understanding of the flow from deeply submerged 


outlets, the following topics are recommended for future study. 
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A detailed investigation on the initial reach (x'/h < 15) 
of the jet. This would be very useful in engineering 
practice, particularly in hydraulic engineering. 
Measurements of the velocity components u, V, Wand the 
turbulent shear stresses for an air wall jet. This would 
check the basic assumptions made in the similarity analysis 
and might give better insight to the mechanics of ‘the flow. 
A satisfactory method of analysis considering the effects 
of the side walls. Study of this kind would in general be 
very useful in practice. 
More experimental observations on swirling wall jets. This 
would supplement the present exploratory studies and cover 


a wider range of swirl numbers. 
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APPENDIX A 
TURBULENT COMPOUND ANNULAR SHEAR LAYERS 


A.1 Introduction 

Consider a circular jet of diameter d issuing from a nozzle 
with an uniform velocity Uo into a large stagnant mass of the same fluid. 
Due to the velocity discontinuity on the periphery of the jet, shear 
stresses are set up and a shear layer originates at the plane of the 
nozzle. In most practical cases, these shear layers become turbulent 
very close to the nozzle itself. This shear layer is termed as a 
simple shear layer. If the surrounding fluid is also in motion, then 
the shear layers are referred to as the turbulent compound annual shear 
layers. A schematic representation of such a shear layer is shown in 
Figure A.1, in which U, is the velocity of the surrounding fluid, x 
is the axial distance measured from the nozzle, x is axial distance 
measured from a virtual origin, r is the radial distance, ry and i) 
are the radial distances to the inner and outer edges of the shear 
layer, u is the turbulent mean velocity in the axial direction, 


Us =U 


i 0° U = u-U 


1 and U. = Uy-Uy - 

A similarity type of analysis is presented to show that ry 
and ro vary linearly with x. This greatly simplifies computation in 
the existing relationship derived by Squire and Trouncer (1944). These 


predicted linear variations were then checked experimentally for four 


values of jet velocity to secondary stream velocity ratios. The be- 
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haviour of constant c in the mixing length relationship (Eq. A.6) is 


also discussed. 


A.2 Existing Work 


The simple annular shear layer has been studied by Kuethe 
(1935), Corrsin (1953), Abramovich (1963), Albertson et al (1950) and 
more recently by Sami et al (1967). Similarity analysis as well as 
experimental observations have shown that both ry and ro vary linearly 
in the axial direction. The velocity distribution in the shear layer 
can be represented by a Gaussian expression (Albertson et al), cosine 
expression (Squire and Trouncer) or by the expression developed by 
Schlichting (1968) for the turbulent far wake. 

The compound annular shear layer has received only very 
little attention. On the experimental side, Abramovich quotes the 
results of Yakovlovskiy for A = U/U5 = 0.4, for an isothermal case 
which showed that for this particular value of the velocity ratio 
the velocity distribution was similar and can be represented by the 


expression, 


u-U + 2 
U (ete — 3/2 
—=7——= [1-7 ] (A.1) 
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r-r 
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Squire and Trouncer developed an integral method for pre- 


dicting the growth of the compound annular shear layer. Assuming 
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co (Eq. A.7) to be a constant quantity, results have been presented for 
various values of A. However, the main predictions of Squire and 
Trouncer have not been checked experimentally. The length of potential 
core given by their method does not agree with the experimental re- 
sults of Forstall and Shapiro(1950). 

The expressions derived by Squire and Trouncer are briefly 
indicated below, as these results would be needed in computing ry 
and iw) by the simplified procedure suggested in the present analysis. 
The equations of motion for the axisymmetric jet can be 


written as 


ou ou wi Vien b9. 

UA Neer ast ako. (A.2) 
3 3 43 
3X (ru) ar or (rv) = 0 (A 3) 


where v is the turbulent mean velocity in the radial direction, t 
is the turbulent shear stress and p is the mass density of the fluid. 
The viscous shear stress has been neglected in Eq. A.2. 


Multiplying Eq. A.2 by pr and on integration, 


r Yr 
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It was assumed that the velocity distribution in the annular shear 


is described by a cosine function of the form 
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where 


Prandtl's mixing length was assumed to be proportional to the width 


of the mixing region, so that 
g = (5-17) (A.6) 
and the turbulent shear stress is given by, 


Z 2 
t= ple(rs-r))] A) (A.7) 
where 2 is the mixing length and c is an empirical constant. At 

this point, Squire and Trouncer assumed of to be invariant with the 
velocity ratio U,/U, and assigned @ constant value of 0.0067, which 
they obtained from experiments on a circular jet issuing into a stagnant 
environment. 


Equation A.4 was evaluated for r = r, and r =(r)tr,)/2 and 


2 
the following two equations were derived. 


(A.8) 
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where ro is the radius of the nozzle, and 
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From Eqs. A.8 and A.9, ri/¥o and ro/¥o were evaluated in a rather 
laborious manner and were plotted against oe for various values of 


dX. These results are discussed later. 


A.3 Present Analysis 


A similarity analysis was performed on the equations of motion, 


making the following assumptions, 
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p= (4 = Fn) (A.11) 
m 
a = h(n) (A.12) 
m 
72 = g(n) (A.13) 
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where b is a length scale, defined as the radial length from the inner 
edge of the shear layer to the point where U = UL, and v is the radial 
velocity. Regarding Eq. A.12, one could develop an expression for ve 
from Eq. A.3, which would be rather tedious. From the experience on a 
number of free turbulent shear flows (Anwar and Weller, 1969) it appears 
reasonable to assume the form given by Eq. A.12. 

Introducing the above expressions into Eq. A.2, it can be 
Shown that, 
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In Eq. A.14, the primes on ry and b represent differentiation with 
respect to x and the primes on f, h and g denote differentiation with 
respect to n. For similarity of flow, the coefficients of n functions 

on the right hand side of the above equations will have to be independent 


of x. They could be constants or zero. Thus 


a x? (A.15) 
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If one assumes that 


(A.16) 


Then in the latter part of the shear layer it is possible to have 
Similarity with q = sh) 1, whereas in the earlier part, it does not 
appear feasible to obtain simple values for the exponents in Eq. A.16. 
This places a restriction on the similarity relationships. Further, 
if a cosine distribution (Eq. A.5) is adopted for the velocity, it 
can be seen that b = (rjtr,)/2, and therefore V5 as well varies 
linearly with x. 

Let R, = r/vo , Ro = r/o and X = x/¥. Then Eqs. A.8 


and A.9 can be written as 
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where x has been taken to be the same as x, and my and m, are the un- 
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known slopes, substitution of Eq. A.19 into Eq. A.18, with the boundary 
condition that R, = | and Rs =eleat-X = 0,7 woulory led 
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For a given X, the coefficients in Eq. A.1/ and A.20 are constants and 
the two equations can be solved simultaneously for R, and Ro for given 
values of X. 

It was found that the results obtained by this procedure 
agree closely with the earlier computations of Squire and Trouncer as 


shown in Figure A.2. 


A.4 Experiments 


The primary jet was produced by a nozzle of inside diameter 
of 0.92", receiving its air supply from the compressed air line of the 
laboratory. The velocity of the primary jet was varied from I2Z5goato 
225.2 ft./sec. The nozzle was located in a rectangular air duct 9" 
wide, 15" deep and 13'-0" long, with suitable entrance arrangements 
(Figure A.3) and an exist diffuser, at the end of which a Tey Wh eae 
blower was located. Four experiments were conducted with the velocity 
of the secondary stream U, equal to approximately zero, 55.0 av one 
and 94.1 ft./sec. The value of A = U,/U5 was respectively 0, 0.257, 


0.496 and 0./5. 


The velocity distribution in the axial vertical plane was 
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measured with a total head tube of diameter 0.05", connected to a 
G6ttingen type 655 inclined manometer which could be inclined in the 
range of 1:2 to 1:25. The manometer could read within an accuracy 

of + 0.5 mm of water for the higher slope and the accuracy could be 

of the order of + 0.05mm of water for the lowest slope. The static 
pressure in the secondary stream was measured with a streamlined probe. 
In all four experiments, the velocity of the secondary stream did not 
vary to any significant extent in the zone of measurements. There was 
a maximum increase of 3% for the case with } = 0.75. The measurements 
were continued close to the end of the potential core of the primary 
jet in the first three cases. However, in the fourth case, for x/To > 
15.7, long term fluctuations in the differential manometer reading 
over a wide range was observed. Hence measurements were not continued 
beyond x/Vo = 15.7. A typical velocity distribution for A = 0.496 is 
shown in Figure A.4. It can be seen that the non-uniformity observed 
in the secondary stream at the nozzle plane, disappears for larger 
values of x. The significant experimental results are given in Table 


Rel 


A.5 Analysis of Experimental Results 


The velocity distribution data for the four velocity ratios 
were plotted with U/U., versus t in Figure A.5 along the cosine curve 
of Squire and Trouncer. It is seen that in Figures A.5a to 5c, the 
experimental observations agree reasonably well with the cosine curve. 


Fora tien case 0b \) = 0.75, the experimental observations show large 
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FIGURE A.4 TYPICAL VELOCITY DISTRIBUTION 
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FIGURE A.5 (CONTINUED) 
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scatter and some departure from the cosine curve. In the plots of 
Figure A.5, the outer radius ro Was used in the calculation of the 
total width of the shear layer. It was somewhat difficult to locate 
the outer edge of the shear layer and some personal judgement had to 
be used in this operation. To overcome this difficulty, the half-width 
b was computed and the velocity distributions are replotted in Figure 
A.6a to d with U/U_, versus n =(r-1)/b. For all the four velocity 
ratios, the velocity profiles are seen to be similar. In Figure A.6a 
the solutions of Tollmien and Goertler for the circular jet (Abramo- 
vich, 1963) are shown. It is seen that for n < 1.0, the Goertler curve 
is closer to the data whereas for n > 1.0, the Tollmien curve is closer 
to the experimental points. For more accurate use, a curve is drawn 
through the points of Figure A.6a for } ~ 0. This curve described 
reasonably well the velocity distribution for the other values of x. 
The invasion of the shear layer into the primary jet is shown 
in Figure A.7 with r/vo versus x/Tq along with the data of Sami et al 
(1967), and Corrsin (1953) for 4 = 0. The experimental data of this 
work show that ri/¥o varies linearly with x/To: If the virtual origin 
for the ry line is fixed at the section where the straight lines of 
Figure A.7 meet the r/¥o = 1.0 axis, then, it could be seen that 
virtual origin is located behind the nozzle and its distance varies 
from 0.4 ro FO] O to aboutie.3 ro TOWeATS Ue oe 


If x, is the length of the potential core, that is the axial 


0 
length from the nozzle to the section where ry is zero, the variation 


of Xo/Vo> as found by extrapolation from Figure A.7 is shown plotted 
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FIGURE A.6 SIMILARITY IN VELOCITY DISTRIBUTION 
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in Figure A.8 with the velocity ratio A. It is seen that Xq/ Vo in- 
creases” Trom=|0-at-\-=-0—to-aboutezSlat \ 020.75, 
The empirical equation given by Forstall and Shapiro (1950) 


as 


errr 2a (A.21) 


does not agree well with the present results (Figure A.8). 

To find out the variation of the empirical coefficient of used — 
by Squire and Trouncer, the value of of X9/To given by the analysis of 
Squire and Trouncer was matched with the present observations and 
was determined. The variation of <2 with A’ obtained by this procedure 
is shown in Figure A.9. For A = 0, cf is equal to 0.0042 and increases 
to 0.0096 for ) = 0.75. It should be mentioned that Kuethe (1935) 
found of = 0.00496 for X\ = 0, whereas Squire and Trouncer used the 
value of 0.0067 for all values of 2. 

The growth of the length scale b is shown in Figure A.10 in 
a dimensionless form. As predicted by the similarity analysis, b 
increases linearly with the axial distance, with the virtual origin 
located about one radius behind the nozzle. 

If b = C, x, then the data plotted in Figure A.11 show that 
C. decreases linearly with X for A up to 0.75. This variation could 


] 
be represented by the equation 


C. 2 0rt15 [0.96 - al (A.22) 
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The values of C, predicted by the equation of Abramovich (1963) are 


less than the experimental values (see Figure A.11). 

As mentioned earlier, it is somewhat difficult to determine 
the outer limit of the shear layer precisely. The observations on the 
growth of r. are shown in Figure A.12, where the increase is seen to be 
linear with the axial distance. If on and dd, are the angles in hearers 
made by the inner and outer edges of the shear layer with the axial 
direction, the variation of these angles with A is shown in Figure 
Ari3. (From Figtre A.13, 1t 1s seen that Qty decreases linearly from 
DpGwedb le ORtOe2 lead Cavs OU ay Oman Qo decreases from 8.9° at 4 = 0 


to10 75. at i =60. 75. 


Summar. 

A similarity analysis of the turbulent compound annular shear 
layer has shown that the length characteristics of this layer vary 
linearly with the axial distance. When this linear relation was used, 
jt was found that the earlier results of Squire and Trouncer could be 
obtained in a much simpler manner. Based on the experimental results 
for four values of the velocity ratio X from about zero to 0./5, the 
empirical coefficient = was found to vary from 0.0042 at X = 0 to 
0.0096 for X = 0.75 whereas Squire and Trouncer assumed a constant 
value of 0.0067 for of for all the velocity ratios. The length of the 
potential core was found to increase from 10 Yo at A ~,0 ta.25 ro at 
x = 0.75, where Yo is the radius of the nozzle. Further, some experi- 
mental results have also been obtained regarding the variation of the 


angle of the inner and outer edges of the shear layer with A and the 


growth of the half width. 
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THREE DIMENSIONAL FREE JETS 
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APPENDIX B 
THREE DIMENSIONAL FREE JETS 


B.1 Introduction 

In Chapter II growth laws and maximum velocity decay for 
three dimensional free jets were worked out from theoretical consider- 
ations. From dimensional analysis it was shown that if Y A is taken 
as a characteristic length parameter to non-dimensionalise the stream- 
wise distance x, the velocity decay relation would be the same regard- 
less of the nozzle shape. In an attempt to find this, some of the 
available data of Trentacoste and Sforza (1966) together with the 
data on a free square jet investigated presently are analysed. The 
experiments on the free square jet were conducted in flume 2 and the 
velocity field explored with a total head tube 0.05 inches diameter. 


The significant data are given in Table B-1. 


B.2 Similarity 
Velocity profiles taken along planes passing through the 
axis of the nozzle are plotted in the conventional non-dimensional 


From Figure B.] and B.2, the two dis- 


forms , Ua! Umo vs. n.. and u/u vs. n,- 


y, 
tributions are found to be similar for distances greater than 5 to 6 
times the nozzle height and closely follow Goertler's free circular 


jet profile. During measurements it was noted that the jet was some- 
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VELOCITY AND LENGTH SCALES FOR A FREE JET 
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B5 
what unsteady and changes direction after noticeable riers: making 
it hard to get consistent readings for the velocities. However, at 
any particular section the profiles are taken a number of times and 
the one looking consistent was taken as the correct profile. This 
explains the reason for getting some scatter in the experimental data. 
The unsteadiness was believed to be caused by the large recirculating 
eddies in the tank and therefore this is likely to happen with experi- 


ments conducted in ambients of finite extent. 


B.3 Velocity Decay 


Variation of maximum velocity Uno with distance is given in 
Figure B.3, which indicates a linear relationship with the virtual 


origin lying at the outlet section. The present data is represented 
by 
120 


ia 7.25 [7 (B.1) 


i 


and the constant coefficient on the right hand side of this expression 
is comparable with the value 7.0 quoted by Yevdjevich (1966). Also 
from Figure B.3, the difference looks small. 

The available data on various shapes were then studied in 
Figure B.4, by plotting Uno! Yo vs. x'/Y A. The shape includes a slender 
nozzle having aspect ratio as high as 10 and considering the odd shapes 
involved, the agreement among the different data looks remarkable. 


Further it is observed that use of Y A as a characteristic length 
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B8 
rather than hydraulic radius r' (Subramanya, 1967) gives closer agree- 
ment. A mean curve is drawn through the test points and is represented 


by the following set of empirical relationships. 
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For most free jets the virtual origin lies very close to the nozzle, 
so that x and x' are more or less identical. From the last expression, 
which is valid in the fully developed region, we find Uno x 1-0 as 
predicted in Eq. 2.38. For a three dimensional free jet of arbi trary 
shape use of Eq. B.2 is recommended when more detail and precise 
values are not known. 

It is to be noted that the exponents on the right hand side 
of the second and third expression of Eq. B.2 are same as the ones 
appearing for a wall jet. Comparison of Figures B.4 and 4.9 shows that 


the maximum velocity decay for a free jet is faster than the decay in 


a wall jet under identical outlet conditions. 
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B.4 Length Scale b 
ee eee 


Abramovich (1963) indicates that the data on length scale 


of axisymmetric and plane jets can be described by a common equation 


bB ys) 


where, h is the height of the nozzle and would be the diameter for 

an axisymmetric case. Since the aspect ratios are | and ~ for these 
Shapes and represent extreme cases, it might be expected that the 
same relationship be valid for three dimensional jets as well. 
Subramanya (1967) analysed some of Trentacoste and Sforza's data for 
various shapes and coricluded that the above expression remains valid. 
The results of the present investigation are shown in Figure B.5, and 
are compared with Eq. B.3. The growth of bY is found to be faster 
than given by the equation and the virtual origin is located at a 
distance of 3.3h in the downstream direction. In any case, the vari- 


ation is linear and is in accordance with Eq. 2.38. 


B.5 Length Scale b, 


Transverse length scales in the x-z plane passing through 
the axis of the jet, as well as in a plane at 0.5 times the nozzle 
height above the axis were found out. Both sets of values are plotted 
in Figure B.5 and are seen to be in agreement with the by scale, which 
is expected. There seems to be no data available on growth of b, for 


any other shape. One set of values for a rectangular jet of aspect 
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ratio 10, measured in a plane passing through the axis of the nozzle 
has been given by Trentacoste and Sforza. They observe that in a 
certain initial length b, decreases and then starts growing in the 
far field in a manner similar to bY growth. This peculiarity seems 


to prevail for slender nozzles. 


B.6 Summary 


Velocity decay relationship can be described by single curve 
regardless of the outlet shape, if ” A is used as a characteristic 
length. In the fully developed flow zone, Uno & aot in accordance 
with similarity analysis. The length scale bY was found to grow 
linearly with distance and the velocity profiles in the fully developed 
flow were similar, following closely Goertler's curve for an axisymmetric 


caSe. 
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APPENDIX C 
VELOCITY PLOTS AND EXPERIMENTAL DATA 
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FIGURE C,1 NON-DIMENSIONAL VELOCITY DISTRIBUTION IN THE CENTER-PLANE 
-CIRCULAR NOZZLE 
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FIGURE C.3 NON-DIMENSIONAL VELOCITY DISTRIBUTION IN THE CENTER-PLANE 
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FIGURE C.4 NON-DIMENSIONAL VELOCITY DISTRIBUTION IN THE CENTER-PLANE 


- TRIANGULAR NOZZLE 
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FIGURE C.5 NON-DIMENSIONAL VELOCITY DISTRIBUTION IN THE CENTER-PLANE 
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FIGURE C.6 NON-DIMENSIONAL VELOCITY DISTRIBUTION IN THE CENTER-PLANE 
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FIGURE C.8 NON-DIMENSIONAL VELOCITY DISTRIBUTION IN THE CENTER-PLANE 
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FIGURE C.9 NON-DIMENSIONAL VELOCITY DISTRIBUTION IN THE z DIRECTION 
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FIGURE C.10 NON-DIMENSIONAL VELOCITY DISTRIBUTION IN THE z DIRECTION 
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FIGURE C.11 NON-DIMENSIONAL VELOCITY DISTRIBUTION IN THE z DIRECTION 


- RECTANGULAR NOZZLE (e =1.50) 
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FIGURE C.12 NON-DIMENSIONAL VELOCITY DISTRIBUTION IN THE z DIRECTION 
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FIGURE C.14 NON-DIMENSIONAL VELOCITY DISTRIBUTION IN THE z DIRECTION 
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TABLE C-2 
FLOW PARAMETERS FOR CIRCULAR WALL JETS 
DIFFUSING IN CHANNELS OF FINITE WIDTH 
Nozzle Diameter = 0.375" 


Tailwater Depth = 55d 


Width of x! Uno Tom Width of Xs Uno Tom 
channel inches ft/sec Ilbs/sft channel inches ft/sec Ibs/sft 
aie 15 24.2 1,485 35.Uu (as) 20.6 1.34 
(U)=24 .62 s0) eal 0.894 (U)=24.62 s)0) 13.0 Oa 
UTS iy Fi ee, GEE SES) a 9.62 0.475 

730 9.68 0.50 10.0 Vod0 0-39 
9.45 7545) p0732 T30 5.98 0.234 
Tea0 Gal Ze pO aaso 16.0 4.90 Gros 
52.0 520m nO elo 20.0 4.26 yal / 
18.0 4.64 0.132 24.0 3.69 0.085 
21.0 GrSOPerUe lo 28.0 3.40 0.0675 
24.0 3.065) 02100 oy aa(t) erate: U.Uos 
(AIER® 3296) | 0.097 36.0 3.04 0.045 
30-0 = 3.90 0.086 
330 3.84 0.082 
36.0 Selly WOM 


0s) 
eTat JAW sasinta 08 ‘erat wise 
NTOIW 3TIMTA IO 2.1aMnAHO At inte ane | 
“Al€.0 = rehome td prea v i | 
ba2 = ddged v5dswi tet aad» 
mo” on” "x 10 AsbEW mo? on” : 
st2\edT 252\3 eofont fanngta Ste\edi s92\t? 29 
he. 3.08 2.8 "OLE 28.1 SAS. - at: . 
(v.0 0.8f 0.¢ S@.AS=0) hes.0  F.0F at $8.08, 
avsco 60: Sa. 2.4 (onal 260.0 38.0 ~-Oxde a2 \ 3" 
c£.0 av.t 0.0f 02.0 83,@ ov 
hMES.0. 60. #82 0.8 SE.0 @p.¥ 2.2 | 
get10 oer oval SRO Stem. SE Ty 
Th.0 aS.p 0.08 Saf.0. 08.2 0.21 
a80.0 at 0.88 : - - SET.0° B.A 0, 8f 
a¥90.0 Of.€ 0.88 Eff.0 06.8 | 1S 
e200 ‘OF.t 0.86 ee ee 
a0.0 BOLE =O BE ve0.0 ae =O AS 


360.0 08.6 O08 
$80.0 8.6 0.88 
tt0.0. 
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TABLE C-2 (continued) 


Width of x Uno Tom Width of x! Uno Tom 
channel inches ft/sec Ilbs/sft channel inches ft/sec Ilbs/sft 
Ae 200) Zia ~ 04928 19.5 3.96 0.102 
(Uj=24.62 4.0 14.8 0.766 23.5 3.40 0.075 
ft/sec) 6/0 walt .2e obeG 27.5 die Ons 
9.0 8.28 0.343 31.5 2.70) 0.040 
12.0 6.35 0.231 36.0 2.50) & 0.032 
15.5 4.98 0.151 
) 7.0! 250 22.25 1.75 
19.0 45m 0811 
(Uj=24.48 4.12 14.8 0.986 
23.0 3.58 0.082 Ft/ Sec) 
se 6.5 10.9 0.59 
270 3.06 0.058 
9.75 oe Ona) 
31.0 3.0 0.05 
13.25 5.36 §©0.173 
36.0 2.75 0.039 | 
17225 3.98 0.1025 
5.5" 2.5 19.3 1.02 22.0 3.10 0.035 
(Up=24.48 5.0 12.12 # 07603 26.5 2.54 0.039 
ft/sec) 8.0 8.65 0.36 31.0 2.40 0.0325 
‘eee 6.07.» 0821 35.0 x10 +0.020 
16.0 4.68 0.136 
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FLOW PARAMETERS FOR SWIRLING WALL JETS 


TABLE 


c-4 


Diameter of nozzle 


Tailwater Depth = 28D 


0275" 
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ee eS 


Swirl No. 


0 
(Uj*=21.13 
ft/sec) 


0.141 


*= 
(Uy 21.4 


ft/sec) 


inches 


aly 
6.0 
10.0 


u 
mo 


ft/sec 


le 


22,9c 
17.30 
11.90 
Soho) 
6.50 
4.59 


25.0 
18%9 
sere 
10%2 
7.45 


b 
y 
ips 


x! 


inches 


eas 
Y = 0.139 


Opec 
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TABLE C-4 (continued) 


Swirl No. x! Uma bY ye DS a3. 
inches ft/sec ft. inches 4 = 0.139 0.512 
0.265 0 29.6 2 3.0 0.09 0.052 
(Up*=21.4 1.5 24.3 0.0637 4.5 0.123 0.083 
ft/sec) 3.0 16.4 0.063 6.0 0.150 0.142 
4.5 14.4 0.046 8.25 0.210 0.227 
6.0 12.4 0.045 10.5 0.265 0.305 
8.25 9.35 0.051 13.5 0.342 0.380 
10.5 Re 0.065 16.5 0.430 0.480 
13.5 5.65 0.084 
16.5 4.62 0.099 
19.5 3.90 Oeil 
22.5 3.50 Se 
25.5 3.0 os 


*Average velocity at nozzle section 


~ nottsse afsson 36 Wooley 2 


(b) TABLE C-4 (continued) 
Bice 2 eee ee ee ee 

Swirl No. xy Tom D 
inches lbs/sft TC. 

0 2.0 1.74 -- 
Sha -- 0.037 

4.0 lees -- 

Teg 0.85 -- 
9.0 -- 0.069 

11.0 OF 57.5 -- 
eit 0.360 0.105 

20.0 ORZ3d5 -- 
225 -- 0.180 

25°.0 URTos -- 

Sca0 0.0832 -- 

40.0 0.045 -- 

50.0 0.025 -- 

60.0 0.01 -- 
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TABLE C-4 (continued) 


Swirl No. ie: Tom T 
inches lbs/sft TU. 
0.141 15 1.49 020325 
4.5 0.84 0.094 
hess) 0.465 Ond275 
12.0 0.238 179 
18.0 Ors lay.2 0.256 
24.0 0.063 0.338 
30.0 0.038 -- 
vie 0.023 -- 
0.265 (De BS iveo6 0.048 
6.0 0.715 Omaha 
1055 0.272 O22) 
sO 0.128 0.336 
20.0 0.0721 -- 
27.0 0.0312 -- 
34.0 Orvis -- 
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